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Abstract. 


The problem of finding quasi discrete energy 


solutions of the non-relativistic Schródinger 


equation when using convential methods of approach, is considered. The study of applying the 
method of configurational interaction based on associated Laguerre functions with single orbital 
exponent in the spin-orbitals is theoretically discussed in Parts 1 and 2 of the work. Investi- 
gations previously reported (Abstr. 1993, 7425/1958) have been extended to the (2p?)*P state 
of He. The last part of the work consists of both a theoretical and an experimental analysis of 
the structure of negative atomic ions. The method of configurational interaction and some 
semi-empirical methods are then applied to finding bound quasi discrete states in H , He 


and Li 
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Introduction. 


Numerous papers dealing with the solution of 
the two-electron problem in atomic theory have 
appeared since the pioneer work by Hy lleraas 
on the helium atom (1). The two-electron systems 
have been the focus of investigations because of 
their considerable interest in physics and astro- 
physics. The investigations being mostly consen- 
trated on the ground states and on some of the 
lowest excited states in the purely discrete region 


of the energy spectrum have also been of guiding 
importance for treating systems containing more 
than two electrons. On the other hand, the appli- 
cation of the conventional theory to the mixed 
discrete-continuous energy interval of the spectrum 
has been severely restricted by purely mathematical 2 
difficulties. How to avoid this trouble, is one of © 
the main topics to be discussed in the present in- | 
vestigation on the quasi stationary states of the 
spectrum. 

The thesis falls naturally into three parts. In the 
first, attention is given to a brief outline of the 
conventional theory. In order to keep the thesis 
to a reasonable length it was necessary to 
restrict this part in several aspects. It is thus in- 
tended to view quite briefly the basic concepts 
involved for non-relativistic treatments of many- 
electron systems. The second part is devoted to 
studies carried out on quasi discrete states of two- 
electron atomic systems with supplementary dis- 
cussions. Finally, in the last part, attention is di- 
rected to a special study of negative atomic ions. 
This part, however, is subjected to the restriction 
that emphasis is placed on the three lightest ions 
H-- He= apd Jae; 


Partai: 
Conventional theory. 


Introduction. 


A common set of mathematical techniques, such 
as the conventional vector, operator and matrix 
formalism in quantum mechanics, is a very useful 
tool for treating many-electron systems. In this 
part there will be given a thorough indication of 
basic concepts which are necessary in non-relativ- - 


istic treatments of bound quantum mechanical 
systems. 


¡Chapter 1. Mathematical formalism. 


L1 Introduction. 

| The present chapter is devoted to pure mathe- 
matical concepts in the study of electronic struc- 
lure when introducing a vector space of finite di- 
hensions, with emphasis on eigenvalues and eigen- 
lunctions. Throughout the thesis the following 
|.otations will be used. 


a, small letters in italics denote vectors 

Fas, both small and large letters associated 
with the index denote operators 

A, large letters in italics denote matrices 

(ab) is the scalar product of two vectors 

(AB) is the matrix product 


a is the complex conjugate of a scalas number 


The star superscript * denotes a osea the 
jomplex conjugated property, whereas the supet- 
bcript + and the symbol — on the top indicate 
fhe adjoint and transposed property, respectively. 


1.2 Vectors and bases. Representation by matrices. 
The application of abstract vector spaces has 
been an important mathematical tool for solving 
physical problems. In fact, in most problems the 
vector space is a Hilbert space of infinite dimen- 
sions, but, practically, one is forced to approxi- 
mate the problem by discussing a space of finite 
Himensions M. 
A physical system may be defined by a mathe- 
atical object called a state vector. This state vector 
completely describes the behaviour of the system, 
_e., if the configuration of the system changes, 
the state vector changes correspondingly. Any vec- 
tor in the space may be written 


M 
x = u;x; where x;= (u; x) 
imd: 
and E x?=(|x/? (1.1) 
¿=1 


| 

Lrovided the vectors u; form an orthonormal 
set, i.e. (uu) = ôy. The vector set U; forms 
i basis, or it is convenient to say that the set 
selects the principal axes of the relevant coor- 
dinate system. Hence, every vector of an M— 
dimensional space may be expressed as a linear 
lombination of any M linearly independent vec- 
tors of the space. 

The basis can be changed by unitary trans- 
Kormations. Let us in this connection consider 
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two orthonormal bases u; and v;. Since u; is 
a basis, each v; may be superposed by the u;'s 
M 


so that v; Zub. 
v=1 


Conversely, since V; is a 


basis, each u; may be written U;= È v;ay;. If 
J=1 
now substituting the former into the latter, one 
finds that the matrix product (AB)=1, or 
B--A-! and A— B-! where A and B possess 
the elements 2,; and 5j, respectively. 
Furthermore, to each vector x in the space 
determined by the basis u; corresponds a 
vector x* in the dual space determined by the 


M 
basis ut so that x — X x;u; where the x;'s 
Al 
are the complex conjugate of the x;'s. The sets 
ut and u; are then said to form a biortho- 
normal basis where (uj u)= ôy. If uj =u, 
the basis is said to be orthonormal. The vector 
x* is usually called the adjoint of the vector 
x. These two types of vectors correspond to the 

bra and ket vectors in the Dirac notations. 
The vector X is represented by a matrix con- 
sisting of M rows and one column, whereas the 
adjoint vector x* is represented by a matrix 
obtainable from that of x by interchanging rows 
and columns and taking the complex conjugate 


of each element. 


1.3 Linear operators. Representation by square matrices. 


The vector formalism also requires the intro- 
duction of operators. In general, any complex 
dynamic variable can be represented by an ope- 
rator. Firstly, let us recall the operator concept. 
An operator Fə is an entity which, acting 
upon an arbitrary vector X and its adjoint x^, 
generates new vectors y and yt 


172 


As can be seen, the operator acts as a pre- 
multiplier to the vector X and as a post-multi- 
plier to the adjoint vector x*. Next, we will 
consider particularly linear operators which occur 
mostly in physical problems. The linearity of an 
operator is determined by the equations 
Fos (x +y) = E, x + Fog y and F (ax) =aFopX . 
For the sake of survey we add some defini- 
tions. To every linear operator Fo, corre 
sponds a complex conjugated operator Fo, defined 
by the equation E, X* =y*, and a transposed 


F,x=y or x* Ej—y* 
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operator F,p defined from equation (y Fy x) = 


(x F,, y) where x and y are two. different vectors. 
The adjoint operator F,, is then the complex 


conjugate of the transposed operator, and satisfies 


the equation (y* Ej; x) —(x' Fysy). It can be 
shown that the adjoint of the adjoint of a linear 
operator is the operator itself. The operator is 
called a Hermitian if it is equal to its adjoint, 
e he EE — op Linear operators occurring in 
quantum mechanics must be Hermitian because 
they represent real dynamic variables, such as the 
energy, the momentum, the position etc. 

Linear operators and square matrices obey the 
same fundamental laws, and there is thus a de- 
finite correlation between these two important 
mathematical tools in quantum mechanics. A linear 
operator is determined if, and only if, its effect 
on each vector u; in the basis is known. That 
is, the corresponding matrix F with elements 
F; must be known 


Fo w= È u; Fy for 7=1,2...M (1.3) 


Hence, if the matrix F is known in a given 
coordinate system, the operator F,; is completely 
determined. The adjoint matrix F^ is similarly 
the representative of the operator F7; in the dual 
space. Analogous to the definitions of complex 
conjugated operator, transposed operator and 
Hermitian operator we thus have complex conju- 


* 
gated matrices F, transposed matrices F and 
Hermitian matrices 


1.4 Spectral representation of a linear operator. 


If a vector x when operated on by an operator 
Fop, is transformed into a numerical multiple of 
itself, then x is called an eigenvector of F, 


(1.4) 


The numerical value of y is called the eigen- 
value of F,,. It must be noted that the eigen- 
vector and the eigenvalue belong to each other. 
Considering an orthonormal basis u;, the eigen- 
value equation may be written 


F,U;—y;u,;, 12 1,2,...M (1.5) 
where y; is the eigenvalue belonging to the 
eigenvector u;. The representative of F,, is now 
a diagonal matrix with elements y;9;. The 


eigenvalue spectrum of F., may consist of a 
discrete set (negative values) and a continuous 


F,px=yx where y being a constant 


set (positive values). It is to be borne in mind 
that the eigenvalues of a Hermitian operator 
are all real. 

Dealing with such Hermitian operators 
represented by Hermitian matrices (square matrices 
with elements H;;=H;;), we can diagonalize 
these matrices by unitary transformations (Section 
7). Obviously such a procedure is equivalent to 
solving the eigenvalue equation (1.5). The correct 
choice of matrices is thus the exact equivalent 
of the proper selection of basic vectors (the 
principal axes of the relevant coordinate system) 
of the eigenvalue equation. In the special case 
of a real Hermitian operator the representative 
is a symmetrical matrix with elements H;; — H;;. 


L5 Projection operators and other types of operators. © 


Extremely useful operators in quantum mecha- 
nics are the projection operators. This operator © 
is usually defined by the product (u;u7) so that © 

Pop () = (wut) = PA (A) (1.6) - 

Obviously the operator is a Hermitian. Firstly, | 
let Po, (j) act as a pre-multiplier to the vector © 
x= u;x;, we obtain P,,(j/)x=ujx;, atid; 
secondly, it acts as a post-multiplier to the adjoint | 
x' = xu}, and we obtain x* Pa (j) —x;uj. 
That is, in both cases the operator generates the 
projections of the vectors on the j-axis. If the 
operation is to be repeated the result is the same 
Pop (2) = Pop (4), Pop (4) Pos (y) =0 if jE 4. Furt- g 
hermore, such an operator is also said to be 
idempotent when Pos (j) Pos (k) = Pos (1) 84; i. e., | 
its only eigenvalues are 0 and 1. Other impor- 
tant Operators in quantum mechanics are the 
permutation operators, such as the anti-symmetri- 
zation operator As = (N) > (-1)4 P where 

E 


p is the parity of the permutation operator P 
and the sum is to be taken over all N/ per- 
mutations. The parity operator also plays an. 
important róle in quantum mechanics. 


1.6 Commuting operators and simultaneous eigenvectors, 


If two dynamic variables can simultaneously 
take definite Wes rad linear operators 
commute {F., Gos} = Fos Gop — Gop Fog = 0. Furt- 
hermore, if two operators commute, then all 
their eigenvectors are common to both, and, the. 
Operators themselves leaving the state vector 
unchanged, represent constants of the motion. 


Ve 


*hysically, this means that the corresponding 
|uantities can be measured simultaneously. 

We conclude this section by including the 
ollowing theorem. If two operators @,, and 
?, commute, and a third operator y. can be 
und to commute with only one of the former, 
en @op and fop are independent of each other. 
therwise, if Y, commutes with both, the @,, is 
. function of fop, or conversely. 


..7 Diagonalization of matrices by unitary 
ansformations. 
We consider only Hermitian matrices, and it 
ill be useful first to recall the definition of 
unitary matrix U by the equations U* eU 
r eqivalendy (U* U) = (UU*) e. 
Next, let us consider a Hermitian» operator 
F, which expressed in an orthonérmal basis 
; is 

Fy= Du; Fu; (1.7) 

2,J 

here F;; are the elements of the representative 
ermitian matrix F. The corresponding eigen- 
alue equation of F, in a chosen basis V; is 


(1.8) 


he basis is selected to be orthonormal 
(v;v;)=5,. The eigenvectors v, may be ex- 
panded in terms of the orginal basis u;, V; — 
E uz. Substituting this into equation (1.8), 
we obtain the matrix equation (F U) = (DU) 
where D is a diagonal matrix with elements 
y 84;. Since both vz and u; form orthonormal 
‘bases, the matrix U is unitary so that (U Dy: 
Hence, we have F — (UDU*) or D = (U+ FU). 
"The Hermitian matrix F is thus said to be dia- 
.gonalized by a unitary transformation. Obviously, 
if two Hermitian operators commute, then they 
can be diagonalized by the same unitary trans- 
formation. 


Fop Ve = Yk Ve > YR = eigenvalues 


| 1.8 Operators acting on functions. 

In practice, it is more convenient to let an 
i pperator F,, in a continuous system act on a 
‘state function Y, rather than on a state vector 


| 


x;. The operation usually generates a new func- 


tion Q; so that 
Fog Wi= Pi (1.9) 
Many of the previous definitions concerning 
| linear operators acting on vectors may be used 
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straightforward to functions, such as the linearity 
conditions Fos{W(x) + 9 (x)} Fo (9) + Fop Co, 
Fos aw(x)] — a Fop (x) where x is a continuous 
variable. The transposed operator m3 is now 
defined from the integral expression 
Sh (x) Fos G(x) dx = Í p (x) Fux) dx. 

Applying the previous definitions for vectors 


in the space a= X u;a (Ò, b=2u,0 (4) and 


b—E,a-2Z£u;F(j)s(j) we see that ^ (2— 
x F(t, j)a(s). If we now let F, act only on 


the a (i's and the ¿(7)'s which both are functions 
of a continuous variable x, we have the relation 
b(x) = Fpa (x)= [F(x,x)a(x)dx. The func- 
tion Z(x,x is called the kernel of the integral 


operator Fo. 
Let f be any dynamic variable represented by the 


operator Fop. Then we have that the mean value 
of f in a given state 4p, is defined by the expres- 
sion 


f Sfo rada (1.10) 


This mean value is in quantum mechanics fre- 
quently called the expectation value uc Bee Or 
the operator. 

The solutions y; of the eigenvalue equation 
(1,5) which now may be written 


Fig Wi (x) — iG, 7= 12,....M (isi) 


are called the eigenfunctions. Usually, the ope- 
rator now acts on an eigenfunction as a differen- 


tial operator. 


19 The Dirac delia function. 
The Dirac delta function 8 (x — y) is defined 
by the equation as follows 


F= f? «yO dy 


Clearly 8 (x — y) 0 when x + y and tends to 
infinity when x=y. The function is thus not 
uniquely defined at x = y; it is used only as a 
symbolic function. According to the definition 
(1.12) the 8-function has a meaning only when an 
integration is performed over its argument so that 


(1.13) 


(1.12) 


fs (x—)) dy —1 when x — y 


The function plays an important róle in the theory 
of eigenfunctions belonging to the continuous 
spectra of operators. 
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1.10 Orthogonality of eigen functions. 

It is easy to prove that eigenfunctions of a 
Hermitan operator form an orthogonal set. Firstly, 
we notice that the eigenvalues in such a case are 
all real. Next, let us consider the eigenvalue equa- 
tions 

Fayl) =y Y (9), 

Fop P(Q) =” P(q) where y and v (1.14) 
belong to a non-degenerate discrete set of eigen- 
values. 

Multiplying the equations from left by q* (4) 


and W* (4) respectively, then subtracting and 
integrating, we obtain 


SOQ) Fa Yla) dq — Jw Fo (o) dq = 
O ASP (4) v (q) dq (1:35) 


Hence, in view of the Hermiticity of the operator 
Fop, it can be written (y—v) fo (1) Y (q) dq = 
O which is equivalent to [o* (4) Y (q) dq — 0 when 
y v. In the degenerate case independent eigen- 
functions belonging to the same eigenvalue will 
be orthogonal. 


1.11 Sturm-Liouville eigenvalue problems. 

In fact, every eigenvalue problem in physics can 
be based on a Sturm-Liouville equation of the 
general form 


Lal) +y w(x) q —0 (1.16) 


where Lo, (7) is a Hermitian differential operator 
Lap (9) = (fq) — 84. 

It must be noted that all the quantities f, g and w 
in the one dimensional case are continuous func- 
tions of the independent variable x in a definite 
interval a < x < b, the function w (x) being a 
type of weighting function. Usually, in most pro- 
blems singular end points occur, and frequently 
infinite intervals. In every case some singular 
point boundary conditions must be imposed on 
the eigenfunctions; they belong either to the 
discrete or continuous set. Such a Sturm-Liouville 
eigenvalue problem can be reduced to a varia- 
tional form (Section 2.7). 


1.12 Quadratically integrable functions. 

An important condition must be imposed on a 
function before it can represent a physical state. 
This condition is that the integral f * y dr 
extended over the fundamental domain of all the 


independent variables must converge. In this way 
the constancy of the integral in time is ensured. 
The corresponding eigenfunction y is then said 
to be quadratically integrable in the chosen domain. 
Ususally, this domain comprises the whole space, 
and the integration over all space yields the prob- 
abilty that the electron is somewhere which must 
be unity. Obviously, every quadratically integrable 


function can be normalized so that Fy* TEE EU 
Conversely, if a wave function is normalizable, 
then it also fulfils the quadratical integrability. Ad- 
ditional conditions being imposed on an eigen- 
function, are especially considered in a later sec- 
tion (Section 2.4). 


1.13 The principle of superposition of states. 
This principle postulates that any state of a 
system can be considered to be formed by super- 
position of other states, and conversely. This 
implies that an arbitrary wave function yy describ- 
ing a state, may be expanded in terms of an ortho- 
normal set of basic functions q; called eigenstates 
or component states 
Y= È a: Qi (1.12) 
where the expansion coefficients a;— [wt dt 
indicate the probability of finding the system in 
the eigenstate Q;. Furthermore, we have that 
the squared modulus |a;|? of each coefficient 
determines the probability of finding the corre- 
sponding eigenstate in the state Y. The nor- 


malization condition - applied to Y leads to 
2 aes. 


Chapter 2. Non-relativistic theory of 
two-electron systems without external 
fields. 


2.1 Introduction. 


In this chapter a brief survey will be given on 
the two-electron atomic theory where all rela- 
tivistic effects and external fields are ignored. 
In principle, the Hamiltonian operator or the 
corresponding Schródinger equation could pro- 
vide us with quantitative results of the eigen- 
values, if and only if exact solutions could be 
found. However, the mathematical difficulty of 
solving the two-electron Hamiltonian is so great 
that an exact solution in a closed form seems to 


> 


quite unattainable (2, 3 and 4). The quanti- 
tive computations of eigenvalues are thus con- 
ined to approximate solutions. 


42 Hamiltonian operator and the corresponding 

Schrodinger equation. Commuting operators. 
The central operator in the atomic theory is the 
amiltonian usually denoted by H,,. This operator 
5 linear and Hermitian. In the absence of exter- 
ial fields the Hamiltonian for a two-electron 
stem assuming that he nucleus is fixed in space, 


explicitly given by 


2 
Ha~ X Hop (Y) + e?/r,, where 
v=1 


Hp) = — h? V3/2m + 
AC ER eS Z e?[r, v9 


F 


E.) 


t is to be noted that v ?, is the Laplacian ope- 
ator, the second term expresses the spin-orbit 
otential, m is the reduced mass, r, the nucleus- 
electron distance and r,, the inter-electronic di- 
stance. The spin-spin and spin-other-orbit inter- 
ctions are neglected. In the case where the 
system moves in a constant internal field, the 
Hamiltonian is independent on time(3H,,/8/— 0). 
he operator itself is a constant of the motion 
hich is termed the energy E. Such a Hamiltonian 
us frequently called an energy operator in quan- 
-um mechanics. The corresponding time-depen- 
ent Schródinger equation of the stationary states 
may then be written 


ih OW; [0f — Hog Y= Ej Y; (2.2) 


with the eigenvalues E; (energies of the sta- 
tionary states) and the exact eigenfunctions w;= 
p; X exp {—ih E; t} where 9; being a func- 
Hon of the coordinates only, satisfies the equation 
op Q; = E; Qj. Hence, in this case with 
conserved energy there exist solutions y; being 
periodic in time. Apart from the phase factor the 
tationary states of the Hamiltonian are determined 
by the ;'s. 

|. Futhermore, the spectrum of the Hamiltonian 
consist both of a purely discrete and a purely con- 
inuous part, but also a mixed discrete-continuous 
interval. However, it is worth while to bear in 
ind that operators commuting with Hop, repre- 
sent constants of the motion (Section 1.6). Such 
loperators are very helpful in the solution of the 
Schródinger equation. 
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If we now first neglect the spin-orbit inter- 
action yes 6 (7) ES in "thé Hamiltonian 
for each electron (equ. 2.1), all the com- 
ponents of the total orbital angular momentum 
and the total spin-momentum commute obviously 
with H,,. Clearly, their squared quantities com- 
mute with H,, too, and may be regarded as con- 
stants of the motion by given quantum numbers 
L and S, Such a Hamiltonian operation is called a 
LS-coupling scheme. Accordingly, the same holds 
true for the total angular momentum given by 
the quantum number J. 

Secondly, if assuming that e”/r,, in Hop is 
small compared with the total spin-orbital po- 
tential V,, the individual total angular momenta 
jp=h+% commute with 24, and are. seally 
good constants of the motion (//—coupling). 

Thirdly, if now neither e?/r,, nor V, are negli- 
gible, only the above total angular momentum 
with quantum number J commutes with Hop. 
This LS]M-coupling scheme holds rigorously for 
an actual atom with J and M as really quantum 
numbers. However, the LS-coupling scheme is 
a very good approximation to two-electron atomic 
systems. 


2.3 The eigenvalue spectrum of the Hamiltonian. 


In atomic quantum mechanics the interest 1s 
particularly focused on finding quantitatively the 
eigenvalues of the Hamiltonian, i. e. solving the 
wave equation 77,,9;— E;q;. This means that 
we are confronted with the problem of finding 
the eigenvalues and the eigenfunctions simultane- 
ously. It is possible to do this approximately by 
using the Rayleigh-Ritz variational method (Sec- 
tion 2.7). In this section we shall only deal with 
the eigenvalue spectrum of the two-electron atomic 
Hamiltonian. As is well known, the eigenvalues of 
the one-electron atomic Hamiltonian occur as ne- 
gative numbers in the discrete and as positive 
numbers in the continuous part, and, further, these 
parts are nicely separated in the spectrum. It is 
also obvious that the two-electron spectrum con- 
sists of a purely discrete set (total negative eigen- 
values) and a purely continuous set (total positive 
eigenvalues), but, in the interval of these, lies a 
mixed discrete-continuous part (total negative 
eigenvalues). This energy interval will be thor- 
oughly discussed in Part 2 of the thesis. 
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2.4 Conditions imposed on the eigenfunctions. 
For convenience, we consider first the one- 
electron problem, i. e., one electron moving in the 
Coulomb field of a fixed proton. On that as- 
sumption, the problem may be reduced to the 
hydrogenic three-dimensional case in which the 
solution can be done exactly by separating the 
Schrédinger equation into three Sturm-Liouville 
one-dimensional equations. At is generally known, 
the same procedure can only be used approximately 
to atomic structure involving two or more electrons. 
However, it is obvious that we can use the same 
boundary conditions as those imposed on the 
eigenfunction of a Sturm-Liouville eigenvalue 
problem. These conditions are as follows. 


i) Only normalizable eigenfunction can be used 
(quadratically integrable, Section 1.12). 

ii) The eigenfunction must be analytical, i. e. the 
eigenfunction and its derivatives are required 
to be continuous and single-valued at every 
point of the configuration space. It is em- 
phazised that this condition of continuity is 
maintained even for the singular domains. 


iii) Some boundary conditions must be imposed 
on the eigenfunction which make it bounded 
for infinite potential energy and at great 
distances from the nucleus. These conditions 
make in other words the solution of the eigen- 
value equation physically admissible. Clearly, 
no electrons can penetrate into the singular 
domains 7;=0 or r;,—0 (Coulomb hole) of 
the potential energy. From a physical point of 
view we must therefore require that the eigen- 
function is equal to zero everywhere in these 
domains. 

In the case of purely discrete eigenfunctions 
the first requirement of quadratic integrability (7) 
is also sufficient at the boundaries (Z7). 

On the other hand, in the continuous case both 
condition (17) and (77) are to be imposed on the 
eigenfunction, but, the condition (7) cannot be 
used directly as in the discrete case because the 
integral Su* wat diverges. This difficulty is 
conventionally avoided by introducing eigendiffe- 
rentials as representing the eigenfunctions 


x (x) = J Y (x, k) dk (2.3) 


where x and & are the coordinates and the con- 
tinuous energy parameter, respectively. The nor- 
malization condition (7) can now be imposed on 


the eigendifferentials p(x). In quantum mecha- 
nics there are wave packets representing electrons 
that, if free, are initially well localized, for which 
the condition (7) holds. 


2.5 Completeness of the basic set. Orthonormal set. 


In the conventional quantum theory of atoms 
and molecules, the eigenvalue problems are ap- 
proximately solved by expanding the eigenfunc- 
tion in terms of a complete orthonormal basic set. 
It is very important to use orthonormal sets be- 
cause we then have ensured that the members are 
linearly independent. Furthermore, principally, the 
members are chosen from a complete set of func- 
tions: Of course, the set does not need to be 
complete, and on this point there have recently © 
been several discussions. However, experience has - 
shown the practical importance of a complete set, : 
i. e. such a set converges much quicker to the true 
limit. The definition of the completeness of an 
orthonormal set may be deduced in the following 
way. 

Let us first consider an infinite series Yy (x) = 


M 
Y Cy Pv (x) where M is a very large number. The 
v=0 


mutually orthogonal functions (,(x) being nor- 
malizable (Section 1.12) are spread out over 
the infinite one-dimensional interval — oo < æ 
<X +00, Next we consider an arbitrary function 
(x) being normalizable in the same interval. 
If the limit 


lim {J|ylPde— |o) 


(2.4) 


the functions p, (x) form a complete set. Obvi- 
ously, if any member is removed, the remaining - 
set will be incomplete. However, it can be shown | 
that eigenfunctions obtained by solving a Sturm- 
Liouville eigenvalue problem ususally form a 
complete set. It is supposed that the problem in 
question possesses an infinite number of ascending 
discrete eigenvalues. We also require that the 
problem subject to special boundary conditions at 
the end points, can be reduced to conventional - 
variational form. Clearly, the sequence of ascend- - 
ing eigenvalues is infinite. The limit of the Mth 
eigenvalue whose eigenfunction has M interior | 
nodes in the above defined interval of x when M- 
tends to infinity, then is- 3 


oo 
fist de- 2 |c]*=0 


lim | Ym) = 00 (2,5) 

M—> co 
Conversely, if now for a Sturm-Liouville eigen- 
ralue problem with only purely discrete eigen- 
values, lim(yw)- oo for the Mth eigenfunction 
with M interior nodes, the set of eigenfunctions 
ust ble complete. Accordingly, the definition of 
rompleteness can be extended to non-orthonormal 
gets, but, such an extension is beyond the scope of 
this section. 

If we now pass to an atomic problem which only 
ossesses a finite number of discrete eigenvalues, 
the corresponding eigenfunctions do not form a 
complete set. Let us suppose that the problem has 
n additional continuous eigenvalue spectrum as 
the limit of a discrete set. Then thé set of eigen- 
unctions is complete if both the discrete and the 
ontinuous set are included. Such a#%et, is the 
iscrete-continuous hydrogen-like set. An entirely 
iscrete set being complete is the associated La- 
uerre functions with single orbital exponent 
¡(Part 2. Section 4.3). 


2.6 The expansion of an arbitrary normalizable 
wave function. 


Generally, the completeness of an arbitrary 
normalizable wave function may now be formed 
iin the following way. 

Let us consider a dynamic variable f represented 
by a Hermitian operator Fo, whose spectrum of 
eigenvalues is not wholly discrete, but is purely 
|discrete in one range of values and purely conti- 
/nuous in another range. It is to be pointed out 
‘that the two ranges are nicely separated, i. e. in 
‘other words, there is no overlapping between 
‘them as in the case of the one-electron spectra. 
"The completeness then requires that an arbitrary 
‘wave function being quadratically integrable is to 
be superposed by all the eigenfunctions of the 
discrete and the continuous spectra as well 


(AE AO Qo 


where the sum is taken over the range of the 
discrete spectrum and the integral over the whole 
range of f in the continuous domain C of the 
spectrum. It is convenient to write 


> 


H Y (g) = Eve Pr (g) (2.7) 
where Gy denotes a summation over the discrete 


. 


set together with an integration over the conti- 


$ 
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nuous set. This statement of completeness is in 
fact an equivalent of the principle of the super- 
position of states (Section 1.13). According to 
this principle the squared modulus | ¢, ? of each 
coefficient in the series (2.6) or (2.7) determines 
the probability of the eigenfunction with eigen- 
value v in the state «y (7). Further the wave func- 
tion can be normalized, i. e. fore dg= Oy |a?'= 
1 since the eigenfunctions of the operator Fep 
form an orthonormal set, I 0329729 
if either x or v belong to the discrete set, and 
[oio dg=8(x—v) if both belong to the 
continuous set. In either case we have the relation 


"EMI py Y da =O» Cu f Q» Pa dg (2.8) 


2.7 Approximation by variational estimates. 


As pointed out in Section 2.2 the stationary states 
of the Hamiltonian are determined by the time- 
independent Schrödinger equation 77,, p; — E; P; 
where q; are functions of the coordinates only. 
Mathematically it is impossible to solve the two- 
electron equation exactly in a closed form, but 
this does not mean that one cannot obtain approxi- 
mate solutions with any accuracy desired. The 
existence of the eigenfunctions and eigenvalues 
has been discussed by several authors (5 and 6). 
The strongest tool so far for deriving such ap- 
proximations is the Rayleigh variational principle. 
In that case the best approximation makes the 
expectation value of the Hamiltonian Hi 


f y* Hop wy dr an extreme value, i. e. 8 (<Hop >) 
=0. Usually one assumes that the trial y has a 
specific form containing a series of variational 
parameters, the best values of which are then 
determined by the extreme condition. In general, 
it can be proved that the expectation value of any 
Hermitian operator yields an upper bound for the 
exact lowest eigenvalue of the operator, that is 


Er >= Y E dE, 

where FE, is the lowest eigenvalue (exact). 
Consequently, if various normalizable trial func- 
tions are used, the lowest expectation value pro- 
duced is the best approximation to the exact 
eigenvalue. The best value for this uppet bound 
is found with high precision for the lowest eigen- 
value of the helium Hamiltonian (4), and it lies 
very close to and above the experimental value. 
In cases where no experimental values are avail- 
able, one also desires to obtain lower bounds to 
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the eigenvalues so that the accuracy of the ap- 
proximation concerned may be determined. In 
most cases such a lower bound can be found, but, 
unfortunately, it is not as reliable as the upper one. 

If now only linear trial functions (the varia- 
tional parameters appear linearly) are used, the 
problem of the eigenvalue optimization is greatly 
simplified. A standard method for solving eigen- 
value problems in such a way is the Ritz varia- 
tional method. 

The eigenfunction «p is expressed by a linear 
combination of a complete orthonormal set e; 


Y= La, (2.9) 
4 


The coefficients 4; are then determined by the 
Rayleigh variational principle. This leads to a series 
of homogeneous linear equations in which the 
eigenvalue appears as an unknown parameter 

ZiHG|£ —y3a]a, —0 (2.10) 
where Z(j|£) — f wf Host, dr . 

The estimate of the upper bound for the lowest 
eigenvalue is then given by the lowest root of the 
secular equation 


|H(j|£ — 85 |—0 (2.11) 


Next, the coefficients and the corresponding 
eigenfunction can be found from the equations 
(2.10) and (2.9). Apparently, the advantage of 
the combined Rayleigh-Ritz variational method is 
the simultaneous solution of the eigenvalue and 
the eigenfunction. It is to be noted that the deduc- 
tion here is confined to orthonormal set in which 
the q”, are linearly independent. By introducing 
non-linear parameters such as the scale parameter 
which is very useful by the optimization, the pro- 
blem of finding the eigenvalue minimum is much 
complicated. Concerning the higher roots of the 
secular equation it can be proved that the Mth 
root is an upper bound to the Mth exact eigen- 
value of the Hamiltonian. Further details on this 
point is given in Section 2.10. 


2.8 Approximation by perturbation of eigenvalues. 


The present section confines itself to a brief 
survey on the conventional theory of solutions of 
eigenvalue problems referred to as perturbation 
problems. For the sake of simplicity the eigen- 
values are chosen here to be non-degenerate. As is 


well known, there are no major difficulties in 
extending the theory to degenerate eigenvalues. 
In general such problems are based on the 
idea that they can be considered as perturbations 
of closely resembling unperturbed problems which 
can be solved exactly. Let us suppose that the 
unperturbed problem is solved exactly, i. e., the 
eigenvalues y? and corresponding eigenfunctions 


ps of a Hermitian operator 77, are known 


Pop P= 5% 1V¥=0,1,.-+. (2.12) 


Next, let us consider an operator Fop+ € Gas 
which differs only slightly from F,,. That 1s, the 
perturbation e G,, is small, 1. e., the parameter e 
approaches zero for any chosen Gop. The pertur- 
bation problem is then represented by the eigen- 
value equation to first-order 


{ Fos te Gis} (p + eg!)= 
{yo + evi} (ye + €Q)) 


where the perturbed eigenfunction 


(2.13) 


9; Fee 
is imposed by the same boundary conditions as 
the unperturbed zero-order functions q^. Carry- 
ing out all the multiplications indicated in 
equation (2.13), making use of the equations 
(2.12) and letting e go to zero in the resulting 
equation, we obtain 


{ Fos — 15} P= 15 Gos } PS 


Clearly, this procedure requires that neither pl 


(2.14) 


nor y! becomes infinite as € approaches zero. 
Since the unperturbed eigenvalues y? are non- 


degenerate it can easily be shown that the first- 
order perturbation to the eigenvalue is 


y, [92* Gp pdr — G(v|v) (2.151 

Then the first-order perturbation to the eigen- 
function can be determined from equation (2.14). 
Considering now a second-order perturbation 
problem the first-order q! gives the second-order 
perturbation to the eigenvalue in terms of spectral 


representation of the unperturbed Hermitian ope- 
rator F,, 


;2— z| GGd»)/Qz—») (2.16) 


Thus the ordinary perturbation technique gives 
pl, yP, etc, in terms of the spectral represen- 


tation of F,,. Unfortunately, the resulting sums 


are usually very difficult to evaluate because the 
largest contributions to them come from the con- 
inuum set of the P's. A method based on the 


ariational principle such as that of the configu- 
ational interaction, is therefore to be preferred. 


.9 Approximation by the configurational inter- 
action method. 
At present there is only one main method 
roviding practical approaches for the eigenvalues 
nd eigenfunctions of many-electron systems. In 
he literature it is often referred to as the polyde- 
or method, the method of superposition of con- 
igurations or the method of configurational inter- 
action (7 and 8). It was already in the pioneer 
ork on He considered as a useful: approach to 
the two-electron problem (9), but, contrary_to the 
method of correlated eigenfunctions,#t was soon 
found that the series of angular configurations 
converged rather slowly which, of course, led to 
much discouragement concerning its utility in pre- 
‘cise calculations of eigenvalues. In the last decade 
ithe method has undergone a remarkable renascence 
lin quantum mechanical treatments of atoms and 
;molecules as well (10 and 11). So far detailed 
¡studies of the method in its application to the two- 
‘electron systems have been made (12 and 13). 
"The convergency is much improved by the in- 
troduction of natural spin-orbitals (14). In this 
section we will outline the method quite briefly 
¡and prefer to make it in the Löwdin view (15). 


Let us consider a normalizable function ¢ (x) 
of a single electronic coordinate x — (r, s) where 
ir and s are the space and spin coordinates, res- 
pectively. According to the completeness theorem 
in Section 2.5, the function can be expanded in 
terms of a complete orthonormal set q;(x). 
| If now using the hydrogen-like set, the set 9; 
must include both the discrete and the continuous 
set (Section 2.6). We prefer to use a complete set 
which is entirely discrete 


p(x)= > 4; Q; (x) where 
a= J p (x) 9$ (x) dx (2.17) 


The one-electron functions p; (r, 5) are usually 
called spin-orbitals. Next, let us consider a nor- 
malizable two-electron wave function (x, , x,). 
The completeness theorem then requires that 


| 
| 
| 
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oo 
ete Sy» LLL (x,) 9; (x5) 
2,7 = 


where (2.18) 
ay n f Y (x; gx») Qi (x,) p; (x) dx, dx, 
Similarly, we have for a normalizable N-elec- 
tron wave function 


CO 
Vd x) 27 
LATE el 
PACA TET e Pr (SS Pr (xn) 
where (2.19) 
allison) m ¥(%,,..,%y) 


X (23 (TER Pi. (xy) dx, ...dxw 


If now making U(x, ,..,Xw) anti-symmetric 
to satisfy the Pauli principle, we must apply 


the operator A. (N) 12 in Section 1.5. 
The result is 


$(,,...,Xw)— 
& SUE is 
(2.20) 
LACA optem a Det (Pe, Pe, TUN 
A selection of N of the indices £,, ka, . . ., Ry 


is called a configuration, and, if the additional 
condition that &, < ka < ... < Ay is included, 
the selection 1s called an ordered configuration K. 
Denoting the normalized Slater determinant be- 


longing to the configuration K with Yg (X, , .. , Xy), 
the wave function is expressed by 


, Xy) 

(2.21) 
55 Pet and 
Ar= (NR q (kı E) 


vafe gens. Xy) Z Ax dac "M 
where 


dx = (EA Det (Pe, Pe, LS 


As can be seen we have deduced the basic 
theorem of the method of configurational inter- 
action, namely, that every normalizable anti-sym- 
metric wave function can be expressed as a series 
of Slater determinants built up from a complete 
orthonormal basic set of spin-orbitals or one- 
electron functions. The Slater coefficients Ax are 
functions of the ordered configurations only, i. e. 
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AgzF(h, b... Ey), and the k-space de- 
scribed by all the values of the indices ky may 
be considered as a sub-space. Applying now (y 
as a trial function in the Rayleigh-Ritz variational 
method, the coefficients A, can be determined by 
the equations (2.10) assuming that the eigen- 
value y in advance is found from the secular equa- 
tion 

|HQ|K) - y8s« | —0 (2.22) 
where H(J|K)=J Y) Hop bx dr. with dr = 
Ox TA d dxy. It is to be noted that J and 


K are two different ordered configurations. 
Furthermore, it can be shown that the Slater 


determinants associated with all ordered configu- 
rations form a complete orthonormal set 


[Uk be dr— yy with Ag [dk ydr (2.23) 


The normalization condition takes the form 


[Ie gar [Asp (2.24) 
where the sum has taken over all the ordered 
configurations. 

In the exact solution the number of ordered 
configurations is infinite, and, so far no such 
solution can be obtained in practice for want of 
any analytical method for solving the corresspond- 
ing infinite secular equation. One is therefore 
forced to truncate the basis. The number of in- 
dependent configurations can be considerably re- 
duced by using properties of the eigenfunction 
determined by the total spin $ and angular mo- 
mentum L, and by the parity of the state considered. 
Since the operators Ss, Lj, and the parity 
operator are all Hermitian and commute with 
Hop (LS coupling), the series can be restricted to 
those configurations or sums of configurations 
which have the same S and L values, and the 
same parity, or it may be expressed as being con- 
ventionally done in the configurational interaction - 
terminology that the configurations must have 
the same species. - ausim dene: y 
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the number of elements to be computed in the 
symmetric matrix of the real Hop is given by 
T) {1 cT). Since (57) is quickly in- 
creasing with M, especially for large N, the order 
of the secular equation and the number of com- 
puting H(J|K) also become so enormously large 
that the numerical work required is formidable. 
On the other hand, if any high-speed electronic - 
computer of sufficiently large capacity is available, 
it is possible to make several such truncated 
approaches without too much trouble. The diffi- 
culty in evaluating H(J|K) depends on the choice 
of the basic set. 

Usually, the convergency of this conventional - 
procedure is very slow. Of course it can be helped 
by using an additional scaling procedure giving 
the satisfaction of the virial theorem. In principle, 
if the basic set is chosen complete, this theorem 
is automatically fulfilled and no scaling is neces- 
sary. More rapid convergency is also obtained by 
introducing natural spin-orbitals into the basic 
set (14). These orbitals may be defined by diago- - 
nalizing the first-order density matrix B so that 


B —(UDU*) or D—(U* BU) where U is a 
unitary matrix and D is a diagonal matrix. (2.25) 
3^ 


The first-order density matrix B is, of course, 
Hermitian, and the above unitary transformation 
can be done as shown in Section 1.7, by in- 
troducing a new orthonormal basis {¢,} in terms | 
of the conventional one {p;} so that = 
Y Pi Ui or, equivalently in matrix form ¢ =g U, 
q — ¢U* with the two basic vector sets Y a 
The functions ¢; belonging to the new 
are called the natural spin-orbitals associated 


linally, it can be shown that the sequence (2.27) 
Onverges much quicker than the conventional 
ine (2.21). 


¡.10 The separation theorem. 


The theorem concerning the separation of the 
figenvalues of a matrix has been known in the 
llgebra for many years (16). Hylleraas and 
Undheim (17) and later on MacDonald 
| 18) have emphasized its importance for varia- 
ional problems in quantum mechanics. Recently, 
e problem of finding reliable upper bounds for 
e eigenvalues of higher quantum states has been 
liscussed by Lówdin (19). 


Let us first assume that such higher eigenvalues 
»elonging to certain eigenfunctions exist fer the 
amiltonian and that the eigenvalues are non- 
egenerate. However, the result should probably 
oe unaffected by the presence of degenerate eigen- 
ralues. Next, if we now shall apply the variational 
principle to excited states, we must orthogonalize 
Ihe trial function to the exact eigenfunctions of all 
e lower states. Unfortunately, these exact func- 
dons are usually unknown, and a varational pro- 
cedure along these lines will give uncertain re- 
sults. On the other hand, by the aid of the sepa- 
ration theorem we can proceed with clearer con- 
science by way of the secular equation obtained 
when applying the Rayleigh-Ritz variational me- 
hod (Section 2.7). 


The theorem postulates as follows. If a trun- 
-ated basic set with M members is extended to 
M-+1) members, the roots of the initial secular 
equation will separate the roots of the extended 
secular equation. Making now the basic set com- 
plete as M-o, it follows that the roots will 
approach the exact eigenvalues monotically from 
above. Hence, the mth root of a truncated secular 
»quation is always the upper bound to the mth 


exact eigenvalue so that gu > EY The “Corres: 
bonding eigenfunction Wp obtained from the 


truncated set is automatically orthogonal to all 
che lower approximate eigenfunctions wy payers 
2 ,. Consequently, making now by energy opti- 


nization RU E, , then Wi approaches the exact 
igenfunction, and this optimal approximation 
ill be orthogonal to all the lower exact eigen- 
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Part 2: 
Investigation of the quasi stationary 
states of two-electron atomic systems. 


Introduction. 


So far very little theoretical work has been 
done concerning the overlapping interval of the 
discrete-continuous eigenvalue spectrum of the 
atomic Hamiltonian. Obviously, discrete states 
lying higher than the lowest ionization potential 
of the atom are able to interact strongly with states 
in the adjacent continuum of the singly ionized 
atom. As a result of the conventional theory eigen- 
functions of such states are not associated with 
interaction of discrete configurations and their 
combinations only, but contain certain components 
of normalizable wave packet functions for the 
adjacent continuous states. The conventional per- 
turbation technique is thus much complicated by 
the presence of matrix components of the Hamil- 
tonian connecting the discrete and continuous states 
of equal eigenvalues. Unfortunately, the fact of 
the matter is that this degeneracy always is main- 
tained owing to the overlapping continuum. Cer- 
tainly, it is possible to apply the conventional 
method of configurational interaction by introduc- 
ing a complete discrete-continuous basic set, such 
as the hydrogen-like one-electron funtcions. The 
space part of the continuous spin-orbitals can be 
approximately represented by normalizable wave 
packets. The trouble, however, is the complexity 
encountered with evalution and numerical com- 
putation of the present discrete-continuous ele- 
ments of the Hamiltonian matrix. One is forced 
to proceed along other lines, and the present in- 
vestigation deals with such problems. 

The investigation falls into three chapters. In 
the first, attention is given to the well-known 
autoionization process which frequently occurs for 
quasi discrete states. Furthermore, Chapters Á and 
5 deal with configurational interaction studies of 
imperfectly and exactly quantized quasi discrete 
states, respectively. Each chapter ends with a con- 
clusion. 


Chapter 3. The autoionization problem. 


3.1 Introduction. 

In the early thirties spectroscopic studies of 
radiaton-less transition processes in atomic spectra 
called autoionization provided important evidence 
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of the instability of the quasi discrete states (20). 
The importance of the autoionization process in 
the optical region of atomic spectra was first re- 
cognized by Shenstone and Majorana 
(21). Unfortunately, at the same time very little 
interest on a theoretical basis was shown for this 
peculiat atomic process (22). About a decade 
passed before Ta-You Wu el al. started with 
a series of quantum mechanical computations of 
eigenvalues of some quasi discrete states in helium 
and also with rough estimates of their lifetimes 
(23). Present knowledge of quantum mechanical 
treatment of quasi discrete is rather meagre. 

Recently, interest in the field of radiaton-less 
transitions has been renewed in connection with 
molecular spectra, electron scattering and capture 
problems, collision problems of the second kind, 
a-decay problems etc. It would change completely 
the balance of the present thesis if we should out- 
line each branch of this field. Since the auto- 
ionization in atomic spectra is the basic problem, 
it has therefore seemed worth while to give in 
this chapter a rather detailed survey of it in the 
simple case of two-electron systems. 


3.2 Definition of the energy interval. The discrete 
and continuous degeneracy. 


In the two-electron atomic problem we are con- 
fronted with an overlapping of discrete and conti- 
nuous spectra in the energy interval 

E(A*) « E « 0 where E(A+) 

is the ground state energy of the corresponding 
positive ionized atom. It can also be expressed in 
the following way: The discrete levels are im- 
bedded in the continuum of the ionized atom. The 
diagrams of the energy terms of He and He* are 
given in Figure 3.1. As can be seen, the above- 
mentioned energy interval lies between the singly 
ionization limit a and the doubly inoization limit 
b. To the right in the He diagram is as an example 
shown a discrete level, the (252p)*P level lying 
at the same height as the adjacent state (15&p)!P 
in the continuum joining the series of terms drawn 
to the left. The discrete terms of He* are shown 
quite to the right of the figure. 

For the assemblage of Hamiltonian stationary 
states with energies lying in the defined interval 
we now introduce the concept of quasi stationary 
states. The terms quasi continuous and quasi 
discrete states then should fall naturally, even the 
last category is frequently termed by doubly ex- 
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Figure 3.1. Term diagrams for He and Het. 
Quasi stationary levels lie between the two limits 


a and P. 


cited states in the conventional theory. All the 
quasi discrete states are practically two-fold de- 
generate. This degeneracy is due to the adjacent 
continuous states. We assume, of course, that 
other degeneracies such as the spin degeneracy and 
the Coulomb degeneracy, are removed. Hence, 
there is a disintegration probability of the atomic 
system when it is excited in a quasi discrete state 
with both electrons in discrete levels of the one- 
electron problem. It must be borne in mind that 
there are relatively very strong coupling effects 


between these quasi discrete and continuous states - 


of equal energies. As a result, when a fairly long 
time interval has elapsed, most of the discrete 
states will assume the eigenfunctions of the adja- 
cent continuous states by ejecting one electron 


from the atom and making another to jump . 


to a level of smaller energy. In this way the total 


energy of the system is conserved (radiaton-less — 


transitions). However, the quasi continuous states 


1 


are the most stable in the sense that one electron . 
usually moves in the innermost quantum orbit | 
whereas the another moves in an orbit extending - 
to infinity; i. e. the latter has been flying away ' 


from the atom and will not return in the ab- 


the disintegration probability varies with the 


strength of the coupling effects, and that there are - 


| 
| 
sence of external fields. It is to be emphasized that | 


quasi discrete states for which this 
is zero. 


Such processes occur in atomic spectra as auto- 


probability - 


ionization, in X-ray spectra like the well known | 


Auger effect and in molecular spectra as predis- 
sociation. In this chapter we only deal with the 
autoionization problem. T 


w- 


| 
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13.3 Time-independent perturbations in the LS 

coupling scheme. Configurational interactions. 
| Let us consider a second-order perturbation of 
the time-independent two-electron Schrödinger 
equation 


Hos Wa = Ea Wa (3.1) 


rwhere the w,’s are the eigenfunctions of the 
stationary states without the phase factor and Hop 
is the Hamiltonian H,,— H,,-- Hj, in which 
he spin-orbit interaction term is neglected (LS 
oupling scheme). It is to be noted that the un- 
erturbed Hamiltonian | H7, — H5, (1) + Hop (2) 
ind the perturbation operator Hjjzze*|[r,, (equa- 
Hon 2.1). aa oe 

The problem consists now in apftoximately 
alculating the eigenfunctions Wa from the exact 


igenfunctions yọ of the stationary states of the 
nperturbed system 


Ho up = Ej 13 
where 143 =0, (2, ,4) Xo (22, la), 
Hop (N) Px (nx, lv) = ew Qu (ny ly) 
and Eg = & + €, (3.2) 


£ is well known q (7, 7) are the discrete-conti- 
uous hydrogen-like orbitals. 

According to the ordinary perturbation theory 
¡(Section 2.8), we then have that the energy in 
tthe second-order is given by the infinite series 
ias follows 


Ea= Es -]- H' (a|o) + 


E | H' (a | v) |?/ (Ea — E») (3.3) 


[Here the symbol & denotes as usual a summation 
iof the discrete set together with an integration 
¡over the continuous set of unperturbed eigen- 
‘functions whereas the prime on the symbol indi- 


icates that the term v=a is omitted from the sum- 


‘mation or integration over v. 

| . . oO oO LI 

- Considering — and u; as representing a 
uasi discrete and quasi continuous state of the 
inperturbed system, respectively, we may write 
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u: = (m , li) 9 (m, D), 
E; = — Z? ((m)?+(m)?) Rydbergs 
¿=P (n , A) F5), 


E; = { — Z? (4)? + P) Rydbergs — (3.4) 
where F (k, 1) is any wave packet function whose 
asymptotic form may be expressed in regular 
Coulomb functions corresponding to a free elec- 
tron with azimuthal quantum number / and 
energy equal to k? Rydbergs (k is the positive 
continuous variable parameter). The unperturbed 
energy of a quasi discrete state Ef will be two- 
fold degenerate (E; = E^, when nj =1). It must be 
emphasized at once that this degeneracy cannot 
be removed for the perturbed system because of 


the continuous range of E; on both sides of Ej. 
However, the equation (2.14) representing the 
first-order perturbation will take the form 


{Hos — Ea} £a = (9a — Hop) sa (3:5) 
with « as the quasi discrete state considered. 
We now expand g, in terms of the complete 


orthonormal set of unperturbed functions yy so 


that ga = Ey w uy where G has the same desig- 
nation as in (3.3). Substituting this into equation 
(3.5), replacing Hj,” by Esuy, multiplying by 
43. and integrating over the whole space, we obtain 


H'(a| B) —8« f up? ue. dr + ca (Ez — E) (3.6) 


If the adjacent state in the continuum is represen- 
ted by £, it follows from (3.6) that H'(a|B) —0. 
Hence, as a result, the first-order perturbation 
theory is applicable to a quasi discrete state only 
if the matrix element of the perturbation operator 
H',, which connects the state with the adjacent 
continuum state is actually zero. Further, the con- 
ventional second-order perturbation technique 
(3.3) can in general be applied to a quasi discrete 
state only if, in addition to the above mentioned 
requirement, there does not exist any intermediate 
state y which connects « with 8. Unfortunately, 
for most of the quasi discrete states H'(a|B) is 
different from zero, and it is clear that the 
perturbation theory cannot be used in the conven- 
tional form when treating these states. 
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3.4 The disintegration probability due to 
constant perturbation. 

Let us consider the following problem. We 
know that, at the initial instant /—/, the atomic 
system is in a quasi discrete state, and we require 
to determine the probability that the system at a 
subsequent instant /—/,, will be in the adjacent 
quasi continuous state. For the sake of survey we 
denote the exact eigenfunction of the initial state 
by wp with definite energy Ep. The final state 
has the eigenfunction We4ac. With energies 
ranging from Ey to Ec+ac. From the pertur- 
bation point of view it is evident that the prob- 
ability transition will be greatest with those states 
in the continuum having energies very close to 
E,. Furthermore, the transition takes place be- 
tween stationary states of the Hamiltonian, and, 
clearly, this operator must be time-independent; 
i. e. the perturbation is constant. In every case we 
have to compute a perturbed eigenfunction which 
at the initial instant coincides with a definite 
unperturbed quasi discrete state wp. 

We start with the equation (2.2) in which 
we require that H.z (eff.) = Hi (p,9) + H'op(9) is 
time-independent; ie. H'op will act as a small 
constant perturbation, and only stationary solu- 
tions of H,, do exist. The stationary solutions 
of the unperturbed operator Ho, are y; —u; X 
exp {—¢h E; ;) with Hou; — E; s; . Since the 
4; form a orthonormal set, the 
set Y; is also complete and orthonormal. 

We shall now seek the solution of H,; by 


expanding Y in terms of the U/j's 


y = Ey c, (7) Wy 


complete 


(3.7) 


where the expansion coefficients are functions of 
time and E has the same designation as in (3.3). 


0 
Substituting (3.7) into the equation nz op Y, 


making use of the above defined equations, it 
follows that 


; de. rna 
ih E, q PE Evo Hep Wy (3.8) 
If we now first multiply both sides of equation 


(3.8) on the left by w*, and then integrate, we 
obtain 


dc , 
int =O, H (x | 1) 


X exp (Z^! (E; — Ey) 1| cy (3.9) 


Since the matrix elements H' (x|v) do not depend 
on time, it follows that 


P () = (Gf)! Sy H' (x|v) 


x f «9 (A X exp (ih! (Er — Es) dt (3.10) 


We now assume that the system initially is in 
a definite unperturbed discrete state v and that 


lo) 


all the zero-order coefficients c” are time- 


independent and zero except one; i. e. o =8 y). 
Equation (3.10) then takes the form 


cP (9 = (Gh) H' (x | v) 


t 
X $ exp (ih (Ez — Es) 1) dt (3.11) 


assuming that the time interval is limited by 


b/EL<<+<rT where r is the mean life of the 
quasi discrete state v. In this way we obtain the 
first-order amplitudes at any subsequent time / 


cQ (y) a — (3.12) 


H' (x |v) X exp (t7! (Ez; — E») t] — H' (x|v) 
Es — By 


when x + v. 
The probability of transition to the state x is 
given by the squared modulus 


Ie (P = 4|H' |»)? 
X sin? $ {7 (E; — Es) t /(Ez — E)? (3.13) 


when x = v; i. e. the unperturbed energies of the 
initial and final states are different. As can be 
seen, in the time interval limited by the above 
indicated requirement the probability depends es- - 
sentially on two factors, namely the energy dif- 
ference in the denominator and the matrix ele- 
ment in the numerator, The quantity itself will- 


y5 


uantatively be determined by a balance between 
these two factors. Transitions to states x which 
have unperturbed energies appreciably different 
rom E; are negligible if they do not possess strong 
interaction. matrix elements H'(x|v). 

The first-order probability transition per unit 
pf time from an initial quasi discrete state v with 
efinite energy E, to a final quasi continuous 
tate x with continuous variable energy E, may be 
ritten with the substitution ¢= $ h-!(E,— E,) 


ea OP=ah2|H' (|) PGS) (3.14) 
here G(5,7) — (m+)! =? sin? {oz}. The super- 
cript on the energies indicating the unperturbed 
roperties has been omitted. A 

Let us now consider the behaviour of the 


unction G(S,1). For large ¢ we see that 
lim {G(¢,7)} is zero and infinite when ¢+ 0 
¿> oo 


nd ¿=0, respectively. On the other hand, if 
ntegrating over ¢ from — oo to +œ (the conti- 
uous range of E,) we obtain a finite value since 


MET (3.15) 


In fact the function G(€, 1) for large 1 satisfies 
e conditions of the Dirac §-function (Section 
(1.9). Hence, we can write 

na (P=? n| H’(v|x) 78 (6) (3.16) 


which is equivalent to 
|a? 2 = 20h z| H'(v|x) [8 (E, — E.) (3.17) 


This expression is actually zero except for E, = Ey. 
Uu to now we have ignored the second-order 


perturbation amplitudes cO) which approxim- 
ately may be given by 


o» = (H' (v| p) H (y|*) 
x (A(v|x) — A (|x) ) (G5, — Ey) (3.18) 


1 
| 


where 


A (co | x) = (exp {th (Eo — Ex)} — 1)/(Eo — Ex) 


|. The first-order perturbation gives a result only 
lif the matrix element H'(»|x) connecting the 
¡initial quasi discrete state v and the final quasi 
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continuous state x is different from zero. If this 
matrix is equal to zero, the transition can occur in 
second-order through an intermediate state y for 
which both H'(v|y) and H'(y|x) are different 
from zero. 


3.5 Radiation-less transitions. Selection rules 
of autoionization. 


Firstly, it will be noted that emphasis is placed 
upon the first-order perturbation result (3.17), 
and, secondly, only transitions between quasi 
discrete and quasi continuous states of equal 
energies have to be considered. Clearly, there 
occur only radiation-less transitions for which 
E, — E,, and, as we shall see later, the transitions 
take place from v to x, not the contrary way. Ac- 
cording to (3.17) the disintegration probability 
per atomic unit of time (2.42 x 10—17 sec) may 
then be written 

e(v|x) =27 |H' (v|x) |? (3.19) 
Conventionally, this quantity is called the prob- 
ability of autoionization. As can be seen, it de- 
pends solely on the strength of the coupling ef- 
fects of the states v and x. If there now in ad- 
dition are present strong interactions between the 
initial state v and some lower purely discrete states, 
and e(v|y) is the total probability of these tran- 
sitions, the mean life of the quasi discrete state v 
is given by r=[e(v|x) +o(v|y) sec. It is also 
worth while to note that the half width of the 
corresponding spectral line is given by (25) 


Urs. 5919 LOS ^T eem (3.20) 


Let us see how the autoionization process 
works in the case H’(v|~) 4-0. Consequently, when 
a fairly long time interval + has elapsed, a radia- 
tion-less quantum jump takes place from the 
quasi discrete state to a specific quasi continuous 
state lying at the same height as shown by the 
horizontal arrow in Figure 3.1. The selection of 
this state is determined from the behaviour of 
H'(vlx). A contrary transition is impossible on 
the assumption that no external fields are present 
since one of the electrons has already left the 
atomic system. On the other hand, the other elec- 
tron jumps to lower levels (the ground level is 
usually in question if such a quantum jump is 
allowed) because of producing energy necessary 
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to eject the first electron from the system and 
give it the positive kinetic energy k2? Rydbergs. It 
should be borne in mind that the total energy of 
the system is conserved. The requirement that 
one or more of the remaining electrons must 
jump to energetically lower levels is common for 
all types of radiaton-less decomposition processes. 

The spectroscopic evidence associated with the 
autoionization process is in general a broadening 
of lines whose initial quasi discrete levels are 
subject to its effect. This fact and other points in 
the interpretation of such lines will be discussed 
in Section 3.7. Theoretically, this spectroscopic 
effect is explained by equation (3.20). Assuming 
that the quasi discrete level v is subject to auto- 
ionization but is also capable of strong radiative 
interaction with a lower singly excited discrete 
level y so that £(v[X) >> o(v|y)> 0, we see that 
the width of the corresponding line is much in- 
fluenced by the autoionization probability «(v|x) 
since 4, is approximately proportional to the sum 
e(v|x) +o(v|y). 

The selection rules may easily be deduced from 
the behaviour of H'(v|x). It must be emphasized 
that the deduction is restricted to the LS coupling 
scheme only. Since the operators — L25, Sos, Jos 
and the parity operator all commute with the 
perturbation operator H’,,, H'op will have no 
matrix elements connecting states with different 
quantum numbers L, $ and J and with unequal 
parities. In view of this and previous considera- 
tions in this section, we immediately obtain the 
selection rules as follows: 

Ignoring the existence of intermediate states 
(the second-order perturbation) the autoionization 
of a quasi discrete tate can only occur if there 
exists a specific quasi continuous state at the same 
height (one or more electron jumps to lower 
configurations) with 


i) the same parity. The parity of a state is even 
or odd according to whether >; /; is even 
or odd where /; is the azimuthal quantum 
number of the jth electron. The sum has to 


be taken over all the electrons. 


i) AL=0 and A $—0. It isto. be borne. in 
mind that this requirement is valid only if 
the spin-orbit and spin-spin interactions are 
neglected, but, however, this effect is negli- 
gible in the case of light atomic systems. 


iii) A J=0 where Les. 


3.6 Imperfectly and exactly quantized quasi 
discrete states. 


Clearly, it may happen that quasi discrete states 
occur for which H'(vlx)=0. In this case radiaton- 
less transitions into the adjacent continuum are 
forbidden (3.19); i.e. for such states definite 
continua do not exist into which the atomic system 
can go. The corresponding transition probability 
of first-order e(v|x) is actually zero., and the 
mean life 7 is entirely determined by o(v|y) and 
may thus be considerably long. In other words, 
such states are capable of emitting and absorbing 
radiation like purely discrete energy states. Con- 
sequently, the energy spectrum associated with 
these states will consist of very sharp lines. Spec- 
troscopic evidence for such lines will be dis- 
cussed in the next section. Obviously, it would 
be hopeless to draw any definite distinction be- 
tween states for which H'(v|»)=0 and states for 
which H'(v|x) is a neligible quantity. However, 
theoretical evidence has been growing that the 
lifetime of most of the quasi discrete states sub- 
jected to autoionization in He is of the order 
1018 — 10—"* sec, whereas for the singly excited 
states the order is 10—7 — 10— sec (25). 

It seems therefore reasonable to divide the 
quasi discrete states into two categories according 
to either H'(v|x) +0 or H'(v|x)-—0 and call 
them imperfectly and exactly quantized states, res- 
pectively. Remarks on their eigenfunctions will 
be given in the last section of this chapter. 


3.7 Spectroscopic evidence of quasi discrete states. 


As already emphasized in the introduction, evi- ; 
dence of autoionzation in the optical region of ato- 
mic spectra has been obtained by spectroscopic mea- - 
surements (20). From the preceding discussions it | 
is apparent that the autoinonization effect gives rise - 
to broadening of the spectral lines. Observations 
show that most of the lines originating from quasi 
discrete levels are considerably broadened although 
in some few cases they are very sharp [Co m p- 
ton and Boyce, Kruger (20)]. It is to be - 
emphasized that Allen (20)] found certain lines — 
by an arc spectrum study of copper, of greater 
width than that caused by normal Doppler effect. 
Theoretically, such a broadening will occur only 
if e(v[x) >> o(v|y) ; i. e., the emission lines origi- - 
nate from quasi discrete states with very strong 
autoionization effects. Undoubtedly, many of such - 


v 


lines are very difficult to observe on account of 
itheir weakness and diffuseness. It is to be hoped 
\that an improved technique in future will be able 
ito hold promise for detailed spectroscopic studies 
iof these lines. Till then very little is known about 
ithe characteristic features of the autoionization 
¡process. Of course, radiative transitions can occur 
Ibetween different quasi discrete states of the im- 
¡perfectly quantized category. The diffuseness of 
both the initial and final energy level probably 
| produces an additional broadening of the spectrum 
[line because ay is now proportional to 
2e(v|~) - eQ |y). 

More discussion on this point is beyond the range 
of this section. 

On the other hand, as mentionéd: above, there 
is spectroscopic evidence for sharp. (narrow) lines 
in the energy spectrum of quasi dis¢fete states. 
Compton and Boyce (1928) and Kruger 
(1930) have observed helium lines in the vacuum 
ultra-violet with wave numbers 279715 (A 357.5), 
312118 (A 320.4), 310061 (A 322.5) and 311346 
cm~! (A 321.2 A). The first two lines were ob- 
served to be very sharp with widths assumed 
to be of the maximal magnitude of 0.05 4 (26). 
This is equal to the width of about 50 cm— in 
the spectral region considered. The equation (3.20) 
then gives a maximal total transition probability 
including both the upper and lower level of 107? 
sec—. If we now assume that this probability 
concerns transitions from quasi discrete levels to 
singly excited levels with probability of the order 
107—105 sec—!, it remains a total probability for 
the upper level of maximum 10% sec—*; i. e., the 
lifetime of the upper quasi discrete level should 
at least be 10—* sec. From these considerations it 
is clear that most of the helium lines do not 
originate from imperfectly quantized quasi dis- 
crete states but from exactly quantized quasi dis- 
crete states such as the series (2pz4)'?L with 
p—sLH;L—12,5 .... (27). The line 


(2p2)?P — (1:2p)?P, and further discussion on 
this point is postponed to Section 5.8. 

So far the discussion has been related to lines 
whose upper quasi discrete levels lie at a fairly 
| great distance from the doubly ionization limit 
| (this limit is denoted by b in Figure 3.1). What 
happens if these levels lie closer to this limit? 
Certainly, the situation concerning the imperfectly 
quantized states will change considerably because 


320.4 is certainly caused by the transition. 
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of the appreciably decrease of the magnitude of 
the corresponding mattrix elements H' (v|x). Con- 
sequently, the autoionization probability appro- 
aches nearer and nearer to zero, and the imper- 
fectly quantized states will assume more and more 
the character of the exactly quantized states. Hence, 
there is just a possibilty that sharp lines in emission 
and absorption can also originate from these levels. 


3.8 Concluding remarks on the eigenfunctions. 


This chapter ends with some concluding remarks 
on eigenfunctions applicable to quasi discrete 
states. It has been seen in Section 3.3 that the 
conventional perturbation technique in first-order 
is applicable to quasi discrete states only if 
H’ (v|x) is actually zero. This requirement is only 
satisfied for the exactly quantized quasi discrete 
states. Obviously, they can be treated by the con- 
ventional mehod of configurational interaction 
(Part 1, Section 2.9). 

On the other hand, the validity of treating the 
imperfectly quantized quasi discrete states by con- 
ventional approximation methods can be discussed. 
Although purely mathematical difficulties prevent 
the deduction of a new method of approaching 
this problem, the time-dependent perturbation 
theory (Section 3.4), however, has accounted for 
leading to a physical understanding of the quasi 
discrete states. We certainly conclude that the 
autoionization of an imperfectly quantized quasi 
discrete state is described by an eigenfunction as 
being a linear combination of a properly anti- 
symmetric normalizable and closed discrete func- 
tion yp which at the initial instant represents the 
discrete state and an anti-symmetric function 
We+ac normalizable in the sense of a wave 
packet describing the final continuous state 


Y (1,2) =A Wot BWe+ac 

where A and B are functions of time. (3.21) 

At the initial instant s=0 where A>>B, 
the expectation value of the energy is essen- 
tially determined from Wp, at a subsequent in- 
stant + from both Wp and We+ac, and, finally 
essentially from YWc+ ac when A << B. 
Ideally, the best expectation value obtained by 
optimization, would take about the same value for 
every time value (degenerate energy values). Ob- 
viously, concerning computation of expectation 
value it does not matter at what time we do it. 
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If we now choose the initial time ¿=0 we ap- 
proximately can use yp only. Clearly, it is quite 
impossible to use the conventional perturbation 
theory (the second-order equations 3.3) which 
expands wp in terms of a complete set of hydrogen- 
like wave functions. The coefficients of the terms 
involving the adjacent quasi continuous states will 
be infinite large. To neglect these important 
terms will probably displace the expectation value 
considerably. It may be possible to avoid this 
difficulty by using a complete basic set for wp 
which is entirely discrete. This view and also other 
aspects have been considered in the next chapter 
in which an extensive numerical investigation of 
imperfectly quantized quasi discrete states has been 
done by using the conventional method of con- 
figurational interaction. Of course, much work 
remains to be done before we can claim to treat 
these states with equal success in quantum mecha- 
nical way as with the exactly quantized states. 
Future work on this problem must rather be 
focused on the Hamiltonian operator in the hope 
of replacing it by a type of projection operator 
that separates the imperfectly quantized states from 
the adjacent continuum. 


Chapter 4. Configurational 
interaction study of imperfectly 
quantized quasi discrete states. 


4.1 Introduction. 


In this chapter we shall deal with a configu- 
rational interaction study of the imperfectly quan- 
tized states. Numerical solutions of their eigen- 
functions, although they by no means represent the 
exact solutions, are of special interest in the mo- 
lecular theory of excited discrete states lying in 
the interval of predissociation. 

The starting point is to expand wp in a sum 
of Slater determinants built up from a complete 
orthonormal basis consisting of entirely discrete 
spin-orbitals. The investigation has been confined 
to the helium atom and partly the negative hydro- 
gen ion, and a thorough study has been carried 
out for two states, namely the important (2s2)15 
and (2p?)'S states. Results of these investigations 
are given in Section 4.9. Also the states (352)1S, 
(3p?)*S and (34?)!S have been investigated in 


detail when considering the 10x 10 matrix pro- 
blem (Sections 4.6, 4.7 and 4.8). The remainder 
of the problem has been discussed in Section 4.11. 
The chapter ends with a conclusion. 

The central problem of this investigation, 
namely, the spectroscopic interpretation of the so- 
lutions, is discussed in Section 4.10. 


4.2 Configurational interactions. Mathematical 
complexity. 


The method of configurational interaction is 
completely outlined in Part 1, Section 2.9. 
Throughout the thesis we have used the conven- 
tional basic set [qj] of spin-orbitals œ (r,s). As 
previously mentioned, if using the hydrogen-like 
set, the basic set must include both the discrete 
and the continuous one-electron functions in order 
to be complete. The procedure is then equivalent 
in effect to using the conventional perturbation 
theory (Section 3.3). The importance of the terms 
will thus mainly depend on the magnitude of the 
corresponding matrix elements H’ (v|x) and the 
energy differences of the unperturbed states 
E = E, — E% . In the problem of treating imper- 
fectly quantized states where H'(v|x) + 0 and 
E°=0, it leads to mathematical difficulties if con- 
figurational terms from the adjacent continuum 
are to be included. Apparently, such terms must 
be considered as dominating in the configuration 
interaction. 


4.3 The associated Laguerre functions with single 
orbital exponent. 


The associated Laguerre functions 


3 (x) = 92 Ta) se «p{—*| 


of order v and with argument x satisfy the Sturm- 
Liouville equation 


(Y +Y x — g (|v, 0)]y=0. (41) 
where g for given values of v and o is a function 
of the independent variable x solely, and the y's 
are the discrete eigenvalues ranging from y to in- 
finity. According to Section 2.5 of Part 1, the 
set of functions y(x) is complete and entirely 
discrete because there is-a minimum eigenvalue 
with a nodeless eigenfunction, and the infinite 
sequence of ascending eigenvalues satisfies the 
condition (2.5). The completeness of the set is - 
maintained whether the functions belong to the - 


| 
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| 
|generalized type of the (2¢+2)-order [y=n+ 
q+1, v=2q9+2, w=2, orthogonal in ordinary 
space] or they belong to the (27+1)-order 


| [y—74-4, v=2q¢+1, o=1, orthogonal with res- 
| pect to a positive weighting function w(x)=x]. 
| Throughout the investigation we have used the 


former because the matrix elements can then be 
found systematically using the ordinary generating 
forms of the functions (28). Introducing now r 
as the independent variable, the (27+2)-order 
functions may be written as expressed in (28, 7) 
where the orbital exponent a keeps constant. 


4.4 Application of a complete discrete basis. 
Discussion. . 


From the foregoing survey it is apparent that 
terms originating from interaction wath the adja- 
cent continuum are extremely importafit for the 
imperfectly quantized quasi discrete states. We 
have seen the mathematical complexity encountered 
in taking into account the interaction when using 
a complete hydrogen-like basis. 

In general, the continuum is an important part 
of the nature of a state, and this importance for 
the ground state of an atomic system has recently 
been emphasized independently by many authors 
(29). The conventional Hylleraas basis for the 
helium ground state (1) is not complete, and 
Kinoshita (30) has introduced an alternative 
series of Hylleraas’ functions which forms a com- 
plete discrete set. Clearly, the latter set must to 
some extent represent the hydrogen-like continuum. 
The best possible choice of basis for treating the 
imperfectly quantized states should therefore be 
a set containing the associated Laguerre functions 


"with a single orbital exponent (Section 4.3). 


In the present investigation on the !S-states we 
have used the eigenfunction 


Y (r, £3) = È c (m, m, g |n) DY (mg |r) 
m, q 
X DN (mg |nr) P? (cos 0) + at E PY (cos 9) 
X E c(m, n, g |n) {PY (mg |r) 
m> n 


Xx DN (ng | gry) + DY (ng | 1973) DY (ng | 17) 
(4.2) 


Elc(mn, gl + Elcin m | P=1 
mag M=N,Q 


with the normalization condition 


m 
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The normalized radial orbitals are represented by 
the appropriate set of associated Laguerre func- 
tions of order (2¢+2) with a single orbital 
exponent 


DY (ag | yr) = (29)? (a+ c + 1) 1} 
x ((@—g —1)!}2 (amr)? LOt22 (2w) 


X exp {— vr) (4.3) 
where a and 4 are the principal and azimuthal 
quantum numbers, respectively, and y is an ovet- 
all scale parameter. For convenience we have 
used the usual spectroscopic orbital notation 15, 
DESO: reek) etc. for this basis. Of 
course, each member of the set can be expanded 
in terms of the complete orthonormal set of 


hydrogen-like functions le c | za) 
n 


PN (ag | nr) = a, (a | 9) p? (ng | Zrin) 


where Lj denotes a summation over the discrete 
set together with an integration over the conti- 
nuous set. The coefficients a, (a|7) may be written 


(4.4) 


a, (a| g) = f P% (ng | Zrin) DY (ag nr) dr. (4.5) 


The discrete hydrogen-like orbitals for which the 
quantum number 77, is zero, are explicitly 


P (ng | Zrin) = — Qn)? QZIn? (n g) 7 
Xx (1 —g—1)!}!? 


x (2Zr/n) Aa (2Zr|n) X exp{—Zr/n} (4.6) 


The integral (4.5) is evaluated in Appendix A. 
In this way we are able to find numerically the co- 
efficients and their squares for all the membets 
belonging to the discrete set. However, the set 
(4.4) is complete; i. e. E |a,(alg)|?=1. If we 
now subtract X|z,(a|g)|? of the discrete set from 
the unity, the remainder should be associated with 
the continuum. In order to distinguish the discrete 
hydrogen-like set (4.6) from our basis (4.3) we 
use the spectroscopic notation [1s], [25], [35], 
ess [2p], ...., etc. for the hydrogen-like or- 
bitals. 

In fact, the results could have been improved 
considerably by the use of an additional non- 
linear splitting parameter. However, the purpose 
of the present investigation is not to obtain excel- 


76 Erling Holøien 


lent values of the energies, but rather to see how 
the configurational interaction method applied to 
the imperfectly quantized quasi discrete states works 
with the basis (4.3). It is to be hoped that the 
acccuracy of the results obtained is quite sufficient 
for this purpose and that some improvements in 
the eigenfunctions would not lead to conclusions 
significantly different from those found here. 


4.5 Configurational interaction study of 
the 1S-states. 

An extensive configurational study of the singu- 
let S-states of He has been carried out using the 
complete discrete basis (4.3). The investigation 
has mostly been focused on the imperfectly quan- 
tized s( 232) 9S vanda (2p7) AS states. (ELO Owen 
31). In this thesis these two papers will be re- 
ferred to as Paper I and Paper Il, respectively. 
Concerning computational details we refer to 
Paper I. In the Appendix of Paper I there are 
also given all the elements accurately computed 
for the 20x 20°S "matrix. 


4.6 Solution of the 10x 10 matrix. Localization 
of the roots. 


It is worth while mentioning that the calcula- 
tions have been carried out in essentially two 
different approximations. In the first, we have 
included systematically all possible configurations 
arising from principal quantum numbers 1, 2 
andas Ere (ise 2925 (27) (20) iS 
(35)?, 2p3p, (3pY and (3d). In the second, we 
have extended this to include all ordered con- 
figurations arising from the principal quantum 
number 4. 

The matrix diagoalization which is equivalent 
to finding the roots of the secular equation (I, 3) 
has been carried out on NUSSE (electronic com- 
puter). The remainder function R“ (E) E — i Og) 
of the 10 x 10 matrix secular equation for adjusted 
values of y has been tabulated over a sufficient 
range of E in order to localize all the ten roots. 


Besides being a function of Æ and 7, Da (E) de- 
pends on the diagonal element (e) used for 
the reduction of the determinantal equation. 
We are especially interested in the roots of 


RO (E) which are idential with the eigen 
values of the matrix. The behaviour of this func- 
tion is illustrated in Figure 4.1 and Figure (II, 2) 


when y=1.0, and Figure 4.2 when »=0.625. As 
can be seen, the appearance of ten distinct roots 
is clearly demonstrated. The picture will be 
changed a little by the variation of the scale factor 
(Figure LL sls): 


Energy (- 12.5 E) (Ryd) 


Figure 4.1. Reg as a function of E. y = 1.0, (e) = (15). 
Roots: (152, 152s, 1535. 
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Energy (- 12.5E) (Ryd) 


Remainder Function [{R(E)x1 o*| 


Figure 4.2. Rj as a function of E. 7) = 0.625, 
(e) = (3p). Roots: (3532, (3p), (34). 


4.7 Optimization of the roots by scaling procedure. 

Using a truncated set containing an over-all 
scale parameter it is likely that the (a4)? roots 
approximate their eigenvalues best. In the 10 x 10 
matrix problem we therefore focus our attention 
on the roots: (15)?, (25), (2p)?, (35)?, (3p)? and 
Ga). 

According to the separation theorem (Section 
2.10) the roots always give upper bounds to the 
exact eigenvalues when M—> oo, Recently, Lö w- 
din (19) has deduced a criterion for the accuracy 


of a truncated set 7 corresponding to the jth 


~~ 


4 
$ 
A 


| 
Dot R^ of the secular equation. It has been 
hown that the upper bound for the mean-square 
ror in the function is proportional to the dif- 
hrence Re —E;.If now, as in our case, the ma- 
rix elements depend upon the scale parameter 
ne roots will be a function of this. Consequently, 


root of the secular equation the R^^s are always 
ipper limits of E; for every value of the scale 
varameter y (M 1s large). The best approximation 
if a truncated set M is, therefore, that associated 
with the optimum value of ». It is also worth 
mphasizing that the advantage of using a scale 
factor as non-linear parameter is that the virial 
neorem will be fulfilled automatically by opti- 
mization (19). Obviously, we have to find the 
ptimum values nopi separately for each root. 

Let us consider the orthogonality.” problem 
essociated with solutions of imperfectly quantized 
juasi discrete states. The result of Lówdin implies 
hat the absolute values of the overlap integrals 
nay be written 


1 
| vt; 97 47 |< {Re — E? 
k 


A TÍ 


s=1 


(4.7) 


t is to be noted that the y;’s are the lower exact 
Heenfunctions where /—0, 1, 2, ..., (& —1). If 
iow D% (truncated set of order M) is the ap- 
oroximate eigenfunction of an imperfectly quan- 
ized state & corresponding to the &th root of the 
ecular equation, the number of the lower exact 
sigenfunctions is infinite. The latter functions 
rorrespond to the entire infinity of (1s ms)'S 
tates including the continuous states. The series 
»f the right-hand side in (4.7) is infinite, but, 
in every case it can be seen that the series is 
'onvergent because the difference in the demoni- 
nator increases with increasing s. The series then 
»robably possesses a finite sum, and, hence, the 
overlap integrals vanish like {Rp — Ez)! for 
R7» E, by energy optimization. 

| This point has also been discussed in Paper I, 
but it must be confessed that some of the conclu- 
sions drawn at this early stage of the investigation 
are wrong. 

| Figure 4.3 and 4.4 give the energy curves 
a(n) separately for the (1s)*, (25% (2p), (35), 
(3p)? and (3d)? roots. The minima and the cor- 
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responding optimum values of y are given in 
Table 4.1. As a comparison we have given the 


Z 
corresponding values of — and (Z— 5/16)/n. 
n 


Table 4.1. Optimum values of energy and 
scale parameter. 
; Optimum 
Optimum (2525/16) 
E energy (Rhc) d e /n 
(15)? —5.7983 20151 (2:00 (33868 
(21)? —1.5473 0.88, |1.00| 0.84, 
Ope iy 51111968 0.843 |1.00| 0.84, 
(35)? — 0.7001 0:62, 10.67 | 0.56, 
(3p)? | —0.5936 056 11067 ls 0.5.65 
(3d)? — 0.3876 0.337 150267. O 


As can be seen, the scale factor determined 
optimumly for the ground state, agrees with the 


hydrogen-like argument Z whereas the optimum 
n 


values of » for all the quasi discrete states are 
sligthly different from the effective nuclear 
charge with screening constant 5/16 in hydrogen- 
like orbital z. 


1 


104 — 
128 — Scalé Constant 


Figure 4.3. E(n) of the roots (1s), (25? and (2p). 
The energies are given in Rydberg units. 
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Figure 4.4. E(n) of the roots (35), (3p)? and (3d). 


The energies are given in Rydberg units. 


The energy seems to be a very smooth function 
of y for all the roots with no signs of any peculiar 
behaviour. 


4.8 The squares of the optimum expansion coef- 
ficients. The effect upon the expectation value 
of the energy. Discussion. 

We now proceed to a discussion of the behaviour 
of the expansion coefficients. According to the 
principle of superposition of states (Section 1.13) 
the squared moduli of these coefficients determine 
the probability of finding the different eigenstates 
(configurations) in the atomic state in question. 
Figures 4.5 and 4.6 are intended to provide suf- 
ficient illustration of the peculiar behaviour of 
these squares with an accuracy up to and including 
two decimal points. Clearly, these vary with y as 
the expansion coefficients themselves. On the right- 
hand side of Figure 4.6 the optimum values are 
plotted, whereas on the left we have given the 
variation of some of the most important squares 
with 7. The configurations denoted by 1 to 10 are: 
(19%, 1525, (29%, (25)*, 1535, 2535, (35)*, 2p3p, 
(3p)? and (3d). 

First of all, even when the energy varies very 
slowly with y in the neighbourhood of the mini- 
mum, the squared moduli of corresponding con- 
figurations change considerably. For the higher 
roots there seems to be deadly competition between 


163 Jo(aa,alg | 


5 10 
Configura tions 


he 


[nd 


Figure 4.5. | c(m, n, q | ope) * for 
(132, (25)? and (2p)? roots. 
{10 X 10 matrix} 
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1 
* 
i 


107 |e(s,5,2] 0| 


Figure 4.6. 


| c(m, n, q | nopr)|? for the 
(392, (35? and (34)* roots. 
(10 X 10 matrix) 


the most important configurations for taking over 
the leadership. Such a competition is noticeable 
between the configurations (2:)?, 1:2: and (2p)? 
in the (25?) 15 solution, the 2:3; and (35)? in the 
(352) *5 solution and (35)? and (34)? in the 
(3p*) *S solution. In the first case, the principal 
configuration (2s)? wins as expected at the opti- 
mum point, in the second, the two terms seem to 
make a compromise, and in the third, the (34)? 


v 


nfiguration is the superior. The latter is un- 
xpected, but it may be explained in the way that 
he (3d)? term is the only one in the truncated 
t which lies higher than the principal term (3p)?. 
Chis discrepancy will probably disappear as the 
uncated set is enlarged with higher configura- 
ions. 

| The effect of different configurations upon the 
ixpectation value of the energy for optimum scal- 
g is graphically illustrated for the three highest 
ots in Figure 4.7. As can be seen, the (3d)? 
onfiguration in the (3p?) 'S solution makes a 
onsiderable lowering of the expectation value 
hich is probably due to angular correlation 
»ffects. 

To make it more clear we have also tabulated 
e energy contributions in Table 4.2. Perhaps the 
ost striking feature of the data giyen in the 
ble is that all configurations lower than the 
rincipal one tend to raise the energy whereas the 
rontrary is true for all the higher configurations. 
his result increases one's respect for correct in- 
erpretation of the solutions to the spectroscopic 
tates associated with the principal configuration. 
his point will be discussed in more detail in Sec- 
ion 4.10. 
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Configurations 


Figure 4.7. Contributions of different configurations to 


the energy for the optimum 10 7X 10 matrix. 


In the table we have also collected data per- 
taining to optimum values of the squares of the 
coefficients with an accuracy up to two decimals. 
On the whole, there seems to be a relation be- 
tween the magnitudes of these squares and cor- 


Table 4.2. 


Computed data concerning the (352), (35?) and (34?) solutions of the optimum 
singulet 10 X 10 S-matrix. All energy entries are in Rydbergs. 


(352)15 (3p2)15 (342)15 

Terms y = 0.609 y = 0.5625 7] = 0.500 

E |c(m,n,q|)P E |c(m,n,q Dl? E l(cm,n,q))P 

(15)? + 0.00020 0.00 + 0.00162 0.00 + 0.00020 0.00 
[Wey + 0.00042 0.00 + 0.00019 0.00 + 0.00073 0.01 
(25)? + 0.00117 0.075 + 0.00076 0.05 + 0.00189 0.04 
(2p)? + 0.00043 0.01 + 0.00957 0.00 + 0.01343 0.09 
1535 + 0.06630 0.00 + 0.00059 0.00 + 0.00001 0.00 
2535 + 0.09084 0.26; + 0.00153 Onil + 0.00054 0.05 
(35)? — 0.83204 0.26; + 0.08984 OPEL + 0.00254 0.03 
2p 5p + 0.04338 0.13; + 0.40947 0.09 + 0.00821 0.14 
(3p)? 72210.06058....1, 0.25 — 0.66877. 0.19 +0.03411 0.15 
(34)? — 0.00174 . 0.00 — 0.43844 0.39 — 0.44923 0.49 
— 0.70042 1.00 — 0.59364 1.00 — 0.38757 , 1.00 


—————————— 
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responding energy contributions. However, there 
are some appreciable exceptions such as the effects 
upon the energy expectation value arising from 
the 1:3; term in the (3s?2)'5 solution and from 
the 2p3p term in the (3p*)'S solution. In fact, it 
is noticeable that great effects can arise from un- 
important terms. 

Further, it is interesting to investigate the effects 
arising from the adjacent continuum. This resul- 
tant effect can be found by analysing our orbitals 
in terms of the complete hydrogen-like set (Sec- 
tion 4.4). Meanwhile, we postpone this analysis 
to the next section in order to apply it to the re- 
sults obtained by the extended 20 x 20 matrix so- 
lution. We only call attention to the predominant 
rôles played by the principal terms (1s)?, (2p)? 
and (3d)? in the configurational interaction as- 
sociated with their respective solutions as shown 
in Figures 4.5 and 4.6. This fact is also expected 
since 15= (15), 2p — {2p} and 3d= (3d) when 
ness: 


4.9 Solution of the 20 x 20 matrix. Detailed study 
of the (152)*S, (252)!S and (2p?)!S solutions. 
Discussion. 

In this section we are concerned with solutions 
of the 20x 20 matrix; i. e. the previous set is 
enlarged to involve also all possible configurations 
arising from the principal quantum number 4. 
Detailed study of this approximation has been 
carried through for the (152)!$, (252)!$ and 
(2p?)'S states. The major portion of the com- 
putations concerning the first two states is given 
in Paper I, whereas Paper II contains those per- 
taining to the latter. A graphical picture of the 
results is given in Figures 4.8 and 49. As can 
be seen, the main features concerning these ex- 
tended solutions are slightly different from those 
of the 10 x 10 matrix. It is obvious that the op- 
timum scale factor varies with the size of the 
truncated set. The data of Table 4.3 do not 
indicate any particular finite limit of the optimum 
scale factor as the set becomes enlarged. 

It is to be noted that the set of order 4 involves 
the configurations: (1s)?, 1525, (25)? and (2p)?. 

In principle, we know that if the complete set 
were to interact the choise of the scale factor 
would be immaterial, but, in practice, when using 
truncated sets the scale factor must be considered 
as being a purely mathematical non-linear para- 


meter to be determined optimumly in ensuring the 
best possible result. Thus there is no reason for 
visualizing the scale factor physically. 


Table 4.3. Variation of the 
scale factor with the size of the truncated set. 


Order (152) 15 (212) (2p?) 1S 
4 19/6 0.893 0.88 
10 2.01 0.88, 0.84, 
20 2.22 0.94 0.81, 


TON 
pue 


10? |s(a,n,a| 9 |? 
$ 


A= from 4x4 matrix 
B= from 10x10 matrix 
C = fros 20x20 matrix 


10? |o(m,n,q|0)1? 


10 
Configurations 


Figure 4.8. |c(m,n,q|nop |? for the (1:2)!S, (2s2)18 and . 
(2p?)18 solutions (20 X 20 matrix}. 


With regard to the results obtained for the 
ground state solution using the optimum 20 x 20 
matrix, we would point out that equal computa- 
tions have been carried out by Shull and Löw- | 
din (14) with essentially the same results. 

It may be of interest to tabulate the course of 
energy contributions for the (2s?)1§ and (2p?)18_ 
solutions as illustrated in Figure 4.9. This has 
been done in Table 4.4. However, it is to be noted 
that the resulting energies given at different stages - 
of approximation are slightly different from those 
obtained optimumly with truncated sets of the 
same order and including similar terms. d 
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As seen from the table, the basic configuration 
is chosen as the first term. Similar sequences with 
the basic configuration as the first term are also 
used in Figure 4.7. The contributions from terms 
(aq|Bp) of different principal quantum numbers 
may be considerably better if additional non-linear 
splitting parameters are introduced. 

It is an interesting point to examine how the 
adjacent continuum is represented in our solutions 
of the quasi discrete states. We consider only the 
(2:2)'S solution. 

The data of Table 4.5 are calculated from the 
formulae explicitly given for a, (alg) in Appendix 
A. The sum of an infinite number of terms in 
the series of squares of the coefficients is found 
by extrapolation of the computed values. It is to 
be hoped that the limits obtained for the discrete 


E contigurasione p E sets in this way are quite accurate to three decimal 

: S. Ae SY : at í places. As seen, the continuum of the hydrogen- 
Eu E rana 20» 20 Ere like set is most represented in the 3s (26 percent) 
{the (232) and (25?) solutions}. and the 4s orbital (78 percent). The 1s and 2; 


Table 4.4. Course of energy contributions. All entries are given in Rydbergs. 


(252) 1S 
The sequence of configurations: 


(25), (125081327, (2p), 1935, 29°35, - G3) a (3p) (Ye. 1a 454 2545 

354s, (45), 2p 4p, 3p 4p, (Ap? 34 4d, (4d) and (4f 

E : —2.087582, —2.078344, —1.422746, —1.519579, 218516971 55169813 
—1.516986, —1.535653, —1:535846, —1.53592, MS. LEO 
_1.580096, — 1.596597, —1.596821, —1.596827, —1.596835, —1.597031, 
—1.597038, —1.597039. 


O M 
(2p*)'S 
The sequence of configurations: 


(2p), (1s)? 152s, (25)?, 1535, 2535, (35)*, 2p 3p, Gp)", G4), 1545, 2545, 35.45, 

(49°, 2p 4p, 3p 4p, (4p)”, 3d ád, (4d)*, and (4f) 

BE 2. —1.225098, —1.210705, —1.193490, — 1.035268, —1.030796, — 1.030731, 
—1.043190, —1.044916, —1.066408, —1.066453, —1.155222, —1. 201045, 
—1.201055, —1.201057, —1.228961, —1.228961, —1 .229016, —1.233701, 

| — 1.235745, —1.237794, 


——M———————————ee—"T 
——————————————————————— 


82 Erling Holeien 


Table 4.5. Contribution of the hydrogen-like set to our orbitals. y= 1, Z= 2. 


A A M ——— 


The sequence of hydrogen-like orbitals is: (15), {2s}, {3s}, {4s}, {5s}, {6s}, {75}, 185,3 
(T0018 Me "185 SERES 


1s 
la, (a|g)|? : 0.70233, 0.25000, 0.01274, 0.00193, 0.00172, 0.00117, 0.00056, 
0.00026, o al A E * 0.00003 A 
Ya, (a!g)|?: 0.70233, 0.95233, 0.96507, 0.96700, 0.96872, 0.96989, | 0.97045, 
0.97081, tà. Lodo i o. B 097107 eee X 
Discrete set: 0.972, Continuous set: 0.028. 
2n 
la, (a|g) |? : 0.23411, 0.75000, 0.00425, 0.00128, 0.00058, 0.00031, 0.00023, 
0.00012, icy: too re PS HA ai 000000 mm Te au 
La, (alg)|?: 0.23411, ^ 0.98411, ^ 0.98836, 0.98964, 0.99022, ^ 0.99053, 0.99076, 
0,99088, PMR ME ROO ca , . 0.99090, Nor Ud 7 MM 
Discrete set: 0.991, Continuous set: 0.009. 
* yt E 
E ELA M Bs 2, 
— by ido (4 191 


la, (a|g)|: pe 0.00000, — 0.57153, — 0.05153, ^ 0.02225, 0.01246, 


0. 


2/25 tse Qa Tis y teensy 0.00000 — ee RA y 
29. (alg) |?: 0.05203, 0.05203, 0.62356, 0.67509, 0.69 
0.72522, TE "UT. 0.79522, 49 erg ee at 


po Ags A 
i | À q 


SU 


| 


rbitals seem to be mostly represented by discrete 
rdrogen-like functions. A straightforward com- 
itation then shows that terms like 1535, 154s, 
535, and 2545 give the coefficient 0.23 of the 
»ntinuous (1555) terms. In this connection it might 
f= mentioned that equally continuous terms re- 
esented in the ground state solution amount to 
coefficient of 0.14 (32). The above estimated 
iue 0.23 seems therefore to be smaller than 
lhat is expected for this quasi discrete state being 
perfectly quantized. It is obvious, however, 
aper I, Table 1) that terms as 1s5s and 157s 
ter with appreciable coefficients in the con- 
guration interaction. If we now assume that both 
s and 7s contain as much as the 4; orbital of the 
drogen-like continuum, the corresponding coef- 
(cient increases to a value of ca 0.33. This fact, 
ogether with the assumption that other .t;ascon- 
Igurations («—8, 9, ...... ) being equally con- 
inuous hydrogen-like in character will appear with 
onsiderable coefficients, shows that we generally 
nay expect appreciably larger effect due to the 
fontinuum in the case of using a complete set. 


40 Interpretation of solutions corresponding to 
certain spectroscopic states. Discussion. 


As long as the basis does not belong to the 
ets of eigenfunctions associated with the unper- 
arbed system (Part 1, Section 2.8), it may be 
irgued that the interpretation of excited solutions 
torresponding to certain spectroscopic states of 
he two-electron atomic system should be viewed 
ith some suspicion. However, according to the 
iriginal philosophy of the method of configura- 
ional interaction the leading term should be the 
¡CF configuration of the state itself associated 
with its particular species. Other terms would be 
represented ideally by carefully determined SCF 
ronfigurations of all other discrete states, and 
;ontinuous too, being allowed in the configuration 
tateraction. Before proceeding it is emphasized 
hat the term associated with the configuration of 
1e state in question is always called the principal 
term even if it is not the leading one. This con- 
tideration should really hold if a hydrogen-like 
»asis is used; i. e. the leading term is identical 
ith the principal one. 

In view of this, the interpretation of solution 
orresponding to the (157)1S, (2p?)*5 and (34?)?S 
tate is obvious because of the hydrogen-like cha- 
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racter of the principal terms when »=Z/n. This 
conclusion is clearly supported by the behaviour 
of the squares of the coefficients graphically shown 
in Figures 4.6 and 4.8. 

On the other hand, the interpretation of solu- 
tions corresponding to the other spectroscopic 
1S-states is less obvious and depends on subjective 
judgement. Let us in this connection discuss the 
(252)'5 solution. The most important coefficients 
in the present 20 x 20 matrix solution are 0.53835 
(254-048 741 A 0.45217 (19289: 
0.31124 (15)? — 0.30553 (1545) and 0.16227 
(1535) (see Paper I, Table 1). It is easy to find 
the discrete hydrogen-like representation of each 
orbital at the optimum point 7=0.9375 from the 
formulae of the overlap integrals a,(alg) expli- 
city given in Appendix A. The results are 
— 0.5484 for 1s, 0.8332 for 2s and zero for the 
orbitals 2p ,35 and 4s. A straightforward calcula- 
tion then shows that the magnitude of the coef- 
ficient associated with the [2s 2s} term in the 
optimum 20 x 20 matrix solution becomes 0.6739. 
This value is smaller than those of the principal 
terms in the (152)!$ and (2p?)!S solution. These 
coefficients are namely 0.9235 and 0.7751, res- 
pectively (see Paper I, Table 2 and Paper Il, 
Table 1). It thus appears that this smaller value 
is essentially due to the unusual large effect arising 
from the (2p)? term in the (25?)!$ configuration 
interaction. Furthermore, it should be borne in 
mind that the (2p)? configuration is almost 
hydrogen-like in character in the optimum solu- 
tion. Obviously, the appreciable importance of the 
(2p)? term is due to the excessive large angular 
correlation effects as being present in the (252)*S 
state. Similar effects have been quite recently 
noticed in the L-shell of the beryllium ground 
state where the total correlation is theoretically 
proved to be almost purely angular (33). 

From these considerations it follows that the 
interpretation of the present (2s?)*S solution 
should be correct. 


4.11 The remainder problem. The completeness 
of tbe basis. Discussion. 


As previously pointed out (Part 1, Section 2.9) 
the exact solution of a quantum state using the 
configurational interaction scheme, contains an en- 
tire infinity of configurations being ordered and 
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Table 4.6. 
Eigenvalues of the solutions. All entries are in Rydbergs. 


AAA AAA ee 


Theoretical 


Solutions Eigenvalues Experimental data tis 
A er te ee E 

He(15")'s —5:8025 —5.8074 a) — 5.8040 i 

He(29).$ —1.5970 | 9 —1.6984 € —1.5421 

Hlp) S —1.2378. 9 — 1.2658, € —1.2728200.1. 2325 


AAA Li A LL d 


having the same species as the principal term. The 
basis then becomes complete, and all troubles con- 
nected with scaling and optimization will dis- 
appear. On the other hand, in practice we are 
confronted with the mathematical complexity of 
solving an infinite matrix numerically. This may 
possibly be done approximately by using any 
iteration procedure; if using a truncated set the 
effect of the remainder set can be estimated by 
the elementary partitioning procedure (19). 


In the present investigation we do not intend 
to make such an estimation. It must be borne in 
mind that the energy minimization by the scaling 
optimization of the truncated sets reduces con- 
siderably the effects of the remainder sets. This 
fact is clearly demonstrated by the apparent con- 
vergence limits of the energy expectation values 
as graphically shown in Figure 4.9. As already 
emphasized, these limits are upper bounds to the 
exact eigenvalues (M=20). The best possible 
total energies attainable by the 20 x20 matrix 
solution are given in Table 4.6 together with some 
experimental data and results obtained by other 
theoretical computations. 


Some comments may be given on the figures 
listed in the table. Restricting possible transitions 
originating from the quasi discrete states to only 
allowed electric dipole radiation involving two 
electron jumps, probable lower bounds of the 
exact eigenvalues can be computed from observed 
helium lines in the vacuum ultra-violet (20). The 
superscripts in the third column of the table there- 
fore denote as follows: % (2s°)'S — (1s 25)!P, 
2 (2p*)!$ — (153p)*P and 9?(25p*)'S — (15 3p)*P 
associated with the helium lines 279715, 312118 
and 311346 cm”, respectively. 

It is to be emphasized that these indentifications 


are somewhat intuitive and must be considered 
with suspicion. Kruger observed, however, the 
above mentioned lines to be of comparable inten- - 
sity, and no appreciable broadening was detected. - 
Thus it is most reasonable that these lines origi- - 
nate from exactly quantized quasi discrete states. 
In any case, we see that our computed eigenvalues 
may be considered as upper bounds. 


Perhaps the most striking feature of the data 
collected in the table is that the six-parameter 
approximation by Wu and Kou (34) applied 
to the ground state of He gives much better re- 
sult than to the (252)! state. The results are de- 
noted by 4) and c) in the table. This is contrary 
to the results obtained by our 20x20 matrix so- - 
lution. This point is also emphasized in Paper I 
without finding any reasonable explanation. If no 
numerical error has crept into their work this dis- 
crepancy can be explained from a logical point of 
view. The eigenfunction which has been used, is 
of the conventional Hylleraas type. As already 
pointed out, these functions are discrete but do 
not form a complete set. Hence, the adjacent con- - 
tinuum which should be of great importance for 
the (2;2)!S state, cannot take part in the con- 
figuration interaction by using such a set. Com- - 
paratively, this continuum is of. minor importance - 
for the ground state. Approximately, the difference 
in energy of about 0.056 Rydbergs probably 
should originate from continuous term in the 
(252)* solution. 


As a comparison it might be mentioned that | 
the resulting energy contribution from the terms 
1545, 2545, 354s and (45)? in our solution amounts 
to 0.061 Rydbergs (see Table 4.4). A straight- 
forward calculation based on data quoted in. 
Table 1 of Paper I and Table 4.5 shows that 68 
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| 
| percent of this contribution comes from continuous 
hydrogen-like [75 ks} terms; i. e. 0.041 Ryd- 


bergs come from such terms in the configura- 
| tional interaction. 


4.12 Conclusion. 


The reasonable results obtained by the present 
configurational interaction study, make it possible 
that the imperfectly quantized quasi discrete states 
can be treated approximately in the conventional 
way when using a complete discrete basis. The 
merits of using such a procedure lie, among others, 
in the three points as follows. Firstly, the impor- 
tant continuous terms of the hydrogen-like set 
are, to some extent, allowed to take part in the 
configuration interaction. Secondly, the,.use,of the 
Rayleigh-Ritz variational method% guarantees 
results being upper bounds for the exact eigen- 
values. Thirdly, by using the scaling procedure 
the important virial theorem is ensured to be ful- 
filled for each solution. Furthermore, the inter- 
pretation of solutions associated with the higher 
roots of the secular equation seems to be plausible. 
It should be pointed out that the problem is ap- 
parently still far from being solved with equal 
confidence as that of the exactly quantized dis- 
crete states. 


On the other hand, the projection operator 
formalism which has been recently rendered a 
simple and practical tool for treating the spin and 
orbital degeneracies (L ó wdin et al. (19, 35)), 
gives some interesting aspects for a more accurate 
treatment of imperfectly quantized quasi discrete 
states in the future. The problem is, of course, 
that admixture of the eigenfunctions of the quasi 
discrete and continuous states is equivalent to an 
admixture of definite linear combinations of dis- 
crete and continuous radial configurations. The 
quasi stationary states degeneracy is probably due 
to radial effects only. A projection operator acting 
on such a mixture will in principle generate the 
projection of it on the quasi discrete or continuous 
eigenstate. Apparently, an auxilliary Hermitian 


| operator of the two-electron radial type for clas- 


| 


sifying the degenerate states, is needed. In this 
way it may be possible to separate the solutions 
of the quasi degenerate states. The mathematical 


difficulties encountered with such a procedure will 


be left to the future. 
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Chapter 5. Configurational 
interaction study of exactly quantized 
quasi discrete states. 


5.1 Introduction. 


In this chapter we shall deal with a similar 
study of the exactly quantized states lying in the 
quasi discrete interval of the helium energy spect- 
rum. We know that the conventional perturbation 
theory can be applied fully to these states because 
there is a non-interacting term with the adjacent 
continuum (Section 3.8). In addition, it must be 
assumed that very small terms occur interacting 
with neighbouring states of the continuum. In fact, 
there seems, just as well, reliable criterion for 
finding eigenfunctions of these states as for exactly 
quantized states lying in the purely discrete region 
of the spectrum. 

However, it can be easily proved by the con- 
ventional selection rules of autoionization that 
only the series (2572) 5 L, g — L, n > L+1, when 
Ea 2.2 0: belong to the exactly quantized 
states in the quasi discrete interval of the two- 
electron atomic spectrum (27). These states escape 
autoionization. Previous assumptions that also the 
series. (mang) "^L. where x: 2,3, 4,5, ., 
m=x+1, =L>x and n>L+1 belong to 
that category, are wrong (36). 

In the present investigation we have carried 
out a configurational interaction study of the 
lowest (2p?)*P level. The result obtained agrees 
fairly well with the observed value. 


5.2 Separated treatment of the state. Configura- 
tional interactions. 


The absence of an interacting term with the 
adjacent continuum makes a separated treatment 
of the quasi discrete state possible. If we now use 
a complete set of discrete-continuous hydrogen- 
like spin-orbitals, it does not mean that there are 
present non-interacting terms associated with the 
continuum. Similar to states lying in the purely 
discrete region of the spectrum, it may occur that 
such terms make considerable contributions to the 
expectation value of the energy. Clearly, using the 
hydrogen-like basis they cannot be neglected. As is 
well known, taking into account such continuous 
terms, the computational work is much compli- 
cated. 7 
On the other hand, a procedure based on a 
complete discrete set of spin-orbitals analogous 
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to that used for the imperfectly quantized states 
is to be preferred. However, it is to be borne in 
mind that there will only be interaction between 
configurations of equal parity. This implies that 
we have to consider configurations of the type 


(3p), (34)?, 2p4p, 3p4p, (4p)?, 344d, (4d)? 
and (4f)?. Since the state is the lowest level of 
its species, the lowest root of the corresponding 
secular equation is an upper bound to the exact 
eigenvalue. 

Since the matrix elements can easily be inter- 
preted in terms of those previously computed 
(Paper I), we have preferred to use associated 
Laguerre functions of the order (244-2) with 
single orbital exponent for the radial part of the 
spin-orbitals. 


5.3 The choice of variables. 

Attempts to reduce the number of independent 
variables for the two-electron atomic system go 
back to the earlier days of quantum mechanics. 
This reduction was first realized for the purely 
S-states (1) and later on for other quantum states 
(37). About ten years ago a re-investigation of 
more general case was made (38). In the present 
investigation we prefer to introduce the angular 
independent variables defined by Hylleraas. 
These angles, geometrically shown in Figure 5.1, 
obey the symmetric equations obtained by the 
spherical triangle 


cos 9 — cos Ü, cos Y, + sin Y, sind, 
X cos (Pa — Pı) 
cos 0, = cos 9, cos 9 +sin Y, sin Ü cos q 
cos Y, — cos Y, cos UF sin 0, sin Ó cosa 
sin Y, sin & = sin Y, sin 
sin Y sin «=sin Y, sin (P, - Pı) 


sin Y sing =sinY, sin(P,— Pı) (5.1) 


As can be seen, the angles a and « are quite 
symmetric when the electrons interchange. Ac- 
tually, we need only take into account four of 
the seven angular variables with the Hamiltonian 
transformation, either 9, 9, 9, and Q, or Y, a, 


V, and Py when necessary. A perfect symmetry 
is similarly conserved in the volume element. In 
addition, we note that Y commutes with any 
of the components of the total spin SAS 
where S = Spr HiS, and that @ commutes 
with L,, the z-component of the the total angular 
momentum L= L,4-L,. Clearly, Y may be treated 
as constant of the motion with respect to the 
operator Sus, x and @ with respect to Ly := 
i1{0/0p, + 9/09,), (atomic unit, =D 


Figure 5.1. Diagram of the independent variables 


5.4 The eigenfunction. 

Since the Hamiltonian commutes with the square 
operators of the total angular momentum ¡he 
and of the total spin $7, associated with the 
eigenvalues L(L+1) and S(S+1), respectively, 
their eigenfunctions can also be eigenfunctions of 
the Hamiltonian (Part 1, Section 1.6). Mean- 
while, it has been possible to construct two-elec- 
tron atomic angular eigenfunctions F(9, %, p|L, 7) 
of L,¿ from those of given separate angular mo- 
mentum Lis, ı =g (g+ 1) where g= L (38). As to 
the problem of determining the eigenfunctions of 
$2,, it is worth while mentioning that the angular - 
functions fall into two types, namely, the sym- 
metric type including terms as cos [y 9] where 
y=1, 2, ...., and the anti-symmetric type with 
only terms as sin [yp] where y—1, 2, ..... 
Concerning our problem with the (2p?)*P state - 
with L=1 and $—1 we now use the normalized - 
anti-symmetric solution of the angular part 


E(D, 91, P |1, g) = (42)! {3 (2g + 1)/2(g + I}? 
X g (81, P) P; (cos 0) (5.2) 4 


^u 


^ 


| where g(0,, q) — sin 9, sin P = sin Y, sin Y and 
P, (cos #) is the associated Legendre function of 
| first order. The function g(0,, p) is anti-sym- 
metric by electron interchange. Obviously, the 
Hamiltonian eigenfunction can be factorized in a 
| purely anti-symmetric space part and a purely sym- 
| metric spin part with L—1 and $— 1. 

Using this result, we may write the trial total 
eigenfunction of the space part as follows 


W (r1 ra) = Le (m, m q |») 
m, g 


X DY (mg | 71) DN (mg | Nra) F(9, a, > 9 | 1, q) 
+272 E F(0,9,,9|1,2) E c(m,n, q | 1) 
g >n 


X {DN (mg | 17) DY (nglara) 
+ DY (ng | nr) DY (mg | "n e 


with the normalization condition (5.3) 
2|c(m,m,[n)P+21c(m,n,9|mP=1 
mn, q mOn, q 


Finally, it should be noted that the normalized 
radial orbitals are represented by the same basic 
set as in the case of the imperfectly quantized 
states (equ. 4.3). In the present chapter we try 
to keep the same notations as in Paper I. 


5.5 Matrix elements of the Laplacian operator. 

It is convenient to write the matrix elements 
of the Laplacian operator in cartesian coordinates 
as follows 


2 
MG|) e — EJ ws Vey dv 


2 
Ej grad, (U* grad, (Y) dV (5.4) 


where dV — dV, dV, = dx, dy, dz, dx dy, dz . 
In ordinary spherical coordinates r,, Y» and 


Q, we have 


MG) (f Gir n) + 17? fG| Pr Oe) 


v= 


+ "na sin”? O, f Gj | Lv Pv)) dv 


i ow;* dw, 
where f (¿| aß) = z T 
and dv=r? r2 sin 0, sin 9, dr, dr, 


X dà, dd, dp, dp, 


(955) 
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However, the problem is now to transform 
M (|j) in the new variables ry, 9, v, 9, and Q. 
In wiev of the ideal connections between the 
new coordinates and the spherical coordinates 
expressed by equations (5.1), this transformation 
may be easily done if we start from the spherical 
form of M(¿|7) in (5.5). The matrix elements 
then become 


2 
M (|j) == J (f Gi] r n) tro f Gj | 08) 


+ (2—») ri? (f Gj | Orv) + sin? 0, f (j| Pr Gr) 
+ (cotg? 9 — 2 cotg cotg 0, cos P + cotg? 0, ) 
x f Gj| PF) — cotg Dy sin p P (j| 0 p) 
+ cos Ph (j| 0,0) — sing sin 9, P (ij | Opr) 
— cotg Ü sin P P (ij | 6 p) 
+ (cotg 0, — cotgó cos P) sin! 9, h (ij | Pr p)) 
—(-»)5'sn?óf(|]op))4dr (5.6) 


where 4 (7 | ef) = f (| aß) + f Gi | Ba) 


The volume element may now be written 
dt =r? r2 sind sin 9, dr, dra 40 dO, dp de, (5.7) 


We note, however, that the 4j;'s entering the ex- 
pression (5.6) are all scaled functions from the 
basic set (5.3) 


Wi= DN G qn) PY | urs) 
x F(0,9,, P| 1,4) (5.8) 


It can be shown that M(;j) finally takes the 
same form as that evaluated for the *S-configura- 
tions (Paper I, 7 and 16). Computation details 
are given in Appendix B. Furthermore, it can be 
shown that both L(;7j) and e (7|) also adopt the 
the same expressions (Paper I, 7, 8 and 16). 


5.6 Matrix elements of the repulsion energy 
operator. Secular equation. 

The matrix elements L/(;j) are more difficult 
to evaluate in the general case. Fortunately, in 
view of the laborious numerical work, the at- 
tempt succeeded in expressing L’(7|j) in terms of 
the previously computed integrals D;,(]j) in 
Paper I, (equ. 7 and 17). In fact, we get the 
same expression of L'(;j) as in Paper I with the 
only difference that C,,(p|é|g) now may be 
written A 
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(2p + DEA F DA (q+1)! (25 — 29)! 5! (DA 


Css (p|&|g) — (pg)? (Cp + 1)(g4- 1)? (25 + 1)! (s — p)! (s = k)! (5 — 2)! 


(ptytii(pte —y — 1)! 


x y (-1)' 


yeso (p—-y—1»(g—r—-0)!6—ey-1!y 


where ø is the smaller of the two integers p—1 
and g—1. 

The derivation of this expression is viewed in 
Appendix B. The resulting secular equation then 
takes the same form as equation (Paper I, 3). 
The numerical values of both M(7jj) and L(|j) 
can be found from Paper I (Appendix) for the 
(mq, nq) configurations in question, whereas the 
values of L'(;j) are quoted numerically in Ap- 
pendix B. 


5.7 Optimization of the lowest root by scaling. 


A straightforward optimization by scaling can 
be made in this case of the lowest secular equation 
root. The situation is quite similar to that of the 
ground state (1s?)1§. We have carried through the 
investigation in a 10x10 matrix approximation 
using only a truncated set with nembers already 
mentioned in Section 5.2 together with an over- 
all scale parameter. Figure 5.2 shows the energy 
curve E(y). The optimum values are found ac- 
curately to be 5,,,— 0.9297 and Emin=— 1.420207 
Rydbergs. The check of the energy minimum by 
the equations (6) in Paper 1 gives better than 
2 parts in 10* for y and 3 parts in 107 for E. 
Clearly, the eigenfunction fulfils the very im- 
portant virial theorem. 


Table 5.1. Results of the 10 X 10 matrix solution of the (2p?)*P state. 
All entries are in Rydbergs. y= 0.9297 


The sequence of configurations: 


(2pP, 2p 3p, (3p)*, (34)*, 2p 4p, 3p 4p, (4p)*, 3d 4d, (4d)? and (4f)? 


E: —1.4015100, —1.4085102, —1.4103014, —1.4155089, —1.4195725, 
—1.4195737, —1.4195848, —1.4199262, —1.4200745, —1.4202074. 


c(m, n, q | Nopt) 


2 |c(m, n, | Nops)|®: 0.94 + 0.04 + 0.00 + 0.01 + 0.01 + 0.00 + 0.00 + 0.00 


TM———M—— M M—MM—M———————————————— i 


Figure 5.2. E() of the lowest root. Contributions to 


5.8 Numerical results. Comparison of the squares 


The numerical work being mostly carried out 
on the NUSSE computer, gives results condensed 
in Table 5.1. 


: +0.97155, —0.19411, —0.05348, —0.08411, + 0.08865, 
— 0.00036, - 0.00286, —0.01868, —0.01177, 


(5.9) 


the energy for the optimum 10 X10 matrix. 


of coefficients. Discussion. 


0.01230. 


+ 0.00 + 0.00 = 1.00 


v 
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Configurations 
Figure 5.3. |c(m,n,q Nogz)|* for the (252)?P. solution. 
> 

Let us discuss below the data listed in the table. 
Firstly, the best expectation value of the spin- 
independent Hamiltonian operator attainable with 
the 10 x 10 matrix approximation is — 1.4202074 
Rydbergs. This value, however, should repre- 
sent an upper bound to the exact eigenvalue. It 
will be interesting to make an estimate of it. As 
mentioned. previously (Section 4.11), the helium 
lines in the vacuum ultra-violet (20) with wave 
numbers 312118 and 311346 cm originate 
most likely from exactly quantized quasi discrete 
states. If we now accept the allowed dipole transi- 

tion in the LS coupling scheme 
(2p?)*P, — (152p)*Po, it seems to be plausible 
that values like —1.4220 and —1.4290 Ryd- 
bergs will be useful for predicting the exact 
eigenvalues of the (2p?)*P, state from the ob- 
served helium lines mentioned above, respectively. 
It is seen that our calculated value represents 
actually an upper bound to these values. The 
energy contributions from different configurations 
are also plotted in Figure 5.2. We clearly see that 
all the terms succeed in lowering the energy. This 
is also expected for the (2p?)*P state solution. 
If, enlarging the truncated set to be almost com- 
plete in the sense that additional non-linear para- 
meters and configurations would affect the cal- 
culated eigenvalue by less than 0.025 ev, the 
| convergence already established (Fig. 5.2) seems 
to be extremely good (assumed exact eigenvalue 
is —1.4220 Rydbergs). Since both the 2p3p 
and 2p4p terms enter the total eigenfunction with 
relatively small coefficients, it seems reasonable 


to suppose that the major portion of the remain- 
ing 0.025 ev is divided among 4, f and higher 
configurations. It is impossible from the present 
investigation to draw any conclusion about the 
D and F limits. On the other hand, if the exact 
eigenvalue is assumed to be —1.4290 Ryd- 
bergs, the difference will be 0.12 ev, and this 
seems to be too large in compensating by contri- 
butions from the remaining configurations men- 
tioned above. 

The third line of the table gives the coef- 
ficients of the total eigenfunction in our basis, 
whereas line 4 gives the squares. 

The latter qantities are shown graphically in 
Figure 5.3. The results indicate that important 
terms in addition to the principal term are 2p3p, 
2p4p, (3p)? and (34)?. It is seen that the prin- 
cipal term (2p)? is decidedly superior to the others 
(94 percent). A similar situation exists in the 
ground state solution (20X20 matrix), the prin- 
cipal term amounts to 85 percent. In this connec- 
tion it must be borne in mind that the superiority 
of the principal terms usually decreases by the 
enlargement of the truncated sets. In fact, compar- 
ing the (2p?)*P results with those of the (152)15 
state, we see the striking similarity between these 
exactly quantized states. In addition, the effect of 
the hydrogen-like continuum seems to be unim- 
portant since terms like 2p3p and 2p4p enter the 
eigenfunction with small coefficients (1 percent). 


5.9 Conclusion. 


The preceding configurational interaction study 
supports the requirement previously indicated in 
Section 3.8, namely, that the conventional ap- 
proximation methods can be fully applied to the 
exactly quantized quasi discrete states. Which 
method is the superior, of course, can be discus- 
sed. We are tempted to prefer the configurational 
interaction method in the Rayleigh-Ritz variational 
scheme in view of its simplicity and use an 
entirely discrete basis being complete. 


Partes: 
Negative atomic ions. 


Introduction. ^ 

Historically, research on negative atomic ions 
can be traced back to Thomson's pioneer work 
on mass spectra in 1913 (39). A survey of ad- 
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vances in experimental and theoretical research 
on negative ions accumulated up to 1950 has 
been given by Massey (40). Since then renewed 
attention has been drawn to free negative ions 
physics in different fields. Thus, in recent years 
there has been an increased appreciable interest 
in atomic and molecular collision processes in- 
volving negative ions (24). Such processes occur 
frequently in the studies of the ionized layers of 
the upper atmosphere, in molecular physics, in 
plasma physics and in a new technique using 
current negative ion sources for voltage doubling 
of accelerators (41). An excellent review contain- 
ing developments in the research on negative ions 
up to date has been given recently by Bra s- 
comb (42). 

Emphasis is placed on atomic ions in this part 
of the thesis. Chapter 6 is devoted to a fairly 
short account of the conventional negative ion 
theory, and Chapter 7 deals with current experi- 
mental research on negative ions. Finally, the three 
last chapters are devoted to investigations carried 
out on the three lightest ions. 


Chapter 6. Theoretical considerations. 


6.1 Introduction. 


The theory of negative ions may be regarded as 
very tentative at the present stage. Clearly, re- 
garding the negative ion structures as systems 
with a number of quantum states analogous to 
those of the neutral atoms and positive ions, the 
conventional quantum theory can be applied. The 
computation and prediction of properties con- 
cerning negative ion structures should therefore 
involve the solution of the corresponding non- 
relativistic Schródinger equation. One additional 
condition has to be raised, namely, that the negative 
ion considered must exist as a permanent struc- 
ture, or the state in question must be bound. 

This chapter deals essentially with such theore- 
tical considerations. The most dominant effects 
responsible for a bound negative ion structure 
are the angular correlation effects which are very 
difficult to compute. 


6.2 The electron affinity. The existence o 

negative ions. 
Let us first assume that the negative ion exists 
in the ground state. The corresponding Schró- 


dinger equation then has a definite solution, 
or we may say that the eigenfunction of the 
Sturm-Liouville eigenvalue problem generates 
a total eigenvalue being lower than that of 
the neutral atom. If so, we say that the atom has 
a positive electron affinity. This quantity usually 
denoted by EA(Z) may be defined theoretically as 
the potential energy difference between the ground 
state of the neutral atom with a free electron at 
rest at infinity and the same state of the negative 
ion. It can similarly be done, of course, for higher 
quantum states. This affinity is simply the attach- 
ment potential necessary to bind the extra electron 
on the atomic structure, or the detachment poten- 
tial necessary to detach the electron from the 
ionic structure. If the negative ion does not exist 
in the quantum state in question, an almost instan- 
taneous disintegration of the jon usually occurs 
through a radiaton-less transition into the adjacent 
continuum of the neutral atom. This may happen 
even if the system is in the ground state. Such a 
transition occurs in the ground state (1s?2s)*S of  — 
He— because this state is unbound. 

Accordingly, the necessary condition for creating 
a bound state is that the corresponding electron 
affinity of the atomic state in question is positive. 


6.3 Quantum mechanical considerations. 


We will now consider the negative ion struc- 
ture from a quantum mechanical point of view. 
It is supposed that the ion exists and its nucleus 
is fixed in space. The spin-less Hamiltonian is 
then usually written (atomic units) as a sum of 
differential operators corresponding to the kinetic 
energies of the electrons together with the elec- 
trostatic potentials of the electron-nucleus attrac- 
tions and the inter-electronic repulsions 

1 2+1 
Hop (Ion) =~ di 2 [isa TE 


ie 


Z+1 
22 p) PS oda 
O 

We can neglect relativistic effects, mass effects 
and spin-orbital effects as long as we are treating 
the lightest ions, or we may treat these effects 
as very small perturbations, but, it is to be borne 
in mind that when treating ions heavier than 
Ne~ all these effects must be taken into account. 

In the first place we consider a ground state | 
problem. We write the Schródinger equation 
of the neutral atom "Y 


~~ 


An Investigation of Quasi Stationary States... 


Hop (Atom) 9 4— EA4WU A 
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where 
| pee 
Hy (Atom) = — = pang We 
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and E, is the total exact ground state eigenvalue 
of the atom. 

By the capture of an extra electron the equation 
may be written 


Hop (Ton) Yr = Er V; (6.3) 


where E, is the total exact ground state eigen- 
value of the ion. Let us now write 74 = 
Z e 
— Y E;(Atom) and £¡=— È Esflony— V. 
gm i=l A 
where Æ; (Atom) and Z;(Ion) are the binding ener- 
gies of the atomic electrons before and after the 
attachment of the extra electron, respectively, and 
V, is the attachment potential of the extra electron. 
Hence, according to the definition of the electron 
affinity (Section 6.2) we deduce 


EA (Z)=— $ (E; (Atom) — E;(Ion)) + V. (6.4) 
z=1 


Usually the atomic structure is more compact 
than that of the negative ions becaus all the elec- 
trons move in an attractive Coulomb field. Accor- 
dingly, the differences 
LE;(Atom) — E;(Ion)} ate positive. For creating a 
stable negative ion structure we thus must require 
that 


2 
V, E (E; (Atom) — E; (Ion)} (6.5) 
am 
]t is seen that the smallness of the alteration 
of the electronic binding energies in the atom by 
the capture of the extra electron, has a decisive 


| influence on the stability of the negative ion struc- 
ture. Investigation has shown that the attachment 


| 


| potential V, is roughly an order of magnitude 


smaller than 1 ev, or in most cases even negative. 
It is therefore apparent that the force of the field 
of the atom exerting on the extra electron, has a 
relatively short range. The only sources of the 
attachment potential seem to come from the screen- 
ing of the nuclear charge by all the electrons and 
the mutual extra electron-atomic electrons inter- 


2l 


actions. It is to be emphasized that the latter effects 
play a dominant róle in forming stable negative 
ion structure. 


6.4 The problem of bound excited states. 


So far we have confined discussion to the nor- 
mal states of negative ions. It seems plausible, 
however, that the probability of forming negative 
ion structure is larger when the ion is excited in 
higher quantum states. This assumption should 
essentially hold true for the quasi discrete states 
of the energy spectrum. Obviously, the energy 
differences (E,(Atom) —E;(lon)) are much 
smaller in equation (6.5). Hence, one would 
expect V, which is a very small quantity, to be 
enough to create stable negative ion structure. The 
states, of course, must belong to the exactly quan- 
tized category and preferably possess the charac- 
ter of metastability. 

In view of this fact it is nevertheless unlikely 
to expect that the negative ions have a number 
of such bound states lying either in the purely 
discrete energy region or in the quasi discrete 
energy interval. If so, the negative ions should 
also possess energy spectra containing a finite 
number of emission or absorption lines. However, 
quantum mechanical computations based on the 
conventional approximation methods indicate re- 
sults for H— which are somewhat contrary to this 
view. Detailed discussion on this point is post- 
poned to Chapters 8, 9, and 10. 


6.5 Approximations by conventional methods. 


To find a quantum state of a negative ion 
bound or not it is feasible to apply the Ray- 
leigh-Ritz variational method (Part 1, Section 
2.7). As already pointed out, this method guaran- 
tees results always being upper bounds to the exact 
eigenvalues when M is large. Lower bounds to 
the exact eigenvalues can also be found. 

Let us denote the upper limits thus computed 
for the negative ion and the corresponding neutral 
atom by E,(lom) and E,(Atom), respectively. 
Furthermore, we write E(Ioz) and E(Atom) for 
the exact eigenvalues. The inequalities E, > E then 
hold true both for the ion and the atom. On the 
other hand, we know that the mean square devia- 
tion from the exact eigenvalues increases with the 
number of atomic electrons. Clearly, we cannot 
draw any definite conclusion on the stability of 
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the state from the computed energy difference 
E¿(Atom) — E, (lon). It may happen that the 
ionic structure is stable even if this difference is 
negative. 

However, the inequality 

E.(lon) — E (Ion) > E¿( Atom) — E (Atom) 
seems to be true for all the atoms with increasing 
eigenvalue deviation error. As a result, the con- 
dition E,(Atom) — E.(lon) > 0 indicating a 
stable ionic structure must be regarded as probably 
inadequate and uncertain. It is seen that an esti- 
mate of the root mean square value of the eigen- 
value for-the ion and the atom as well, is highly 
desirable. Unfortunately, such a procedure is very 
comprehensive for ionic systems with a large num- 
ber of electrons, and the certain condition left for 
finding bound states from computed variational 
upper limits is 


Eesptı (Atom) — E. (lon) > 0 (6.6) 


where Eerptı is the experimental value of the state 
considered in the neutral atom. Unfortunately, the 
experimental values are mostly unknown, in parti- 
cular for the excited states. Of course, we can 
compute the upper bounds and the lower as well 
by using the variational method with either cor- 
related or configurational interaction trial func- 
tions. In addition, the mean square deviation from 
the exact solution f | Wexact — Wapprox |? dT 

may possibly be investigated. In this way we may 
get reliable results of EA(Z), but, the labour of 
the work required, especially in cases of heavier 


ions than Li—, will be formidable even if elec- 
tronic computers were used. Such works now going 
on for L;— will be pointed out in Chapter 10. 

There is no doubt that the procedure of con- 
figurational interaction is superior to that of cor- 
related eigenfunctions. Particularly, the computa- 
tional work of evaluating energy integrals involved 
in the matrix elements is much easier. One is only 
forced to include spin-orbitals of sufficiently high 
angular quantum numbers in order to obtain an- 
gular correlation effects in Ve which are ultimately 
required for a satisfactory result. À combination 
of the procedures, i. e., using a correlated factor, 
seems to be plausible and will probably give 
fairly good results in view of its simplicity. To a 
first approximation the factor could be chosen as 
the linear form 
{1+ era + Bret... d vruy] where ry 
are the mutual distances between the extra elec- 
tron and the N remaining atomic electrons. The 
approach itself is simple and straightforward, but, 
the practical computation of most of the energy 
integrals may be fairly cumbersome. 

To get estimates of the electron affinities of 
heavier atoms one has thus been forced to apply 
semi-empirical methods. 


6.6 Semi-empirical methods. 


In this section we shall attempt to present a 
short outline of the semi-empirical methods based 
on various extrapolation formulas for estimating 
ionization or detachment potentials of negative 


Table 6.1. Atomic electron affinities. All entries are in ev. 
———————MÀM M i ee ee 
TN. Quadratic Bates — Bacher — Johnson — : 
extrapolation  Moisewitsch Goudsmit ^ Rohrlich Spona 
dt 1. i i 5c aC pts SEE e aale ade dE: 
Be 0.42 —0.22 —0.16 0.82 Stable 
CH 1,32 0.93 1:25 1.21 1.13 t 0.15 
N —0.10 + 0.1 —0.27 —0.74 0.54 Probably stable 
Qu 1.20 1.56 104 1.47 1.465 + 0.005 
Fun (88:44 3.67 3.39 3.62 3.62 +0.09 © 
AD 1.01 0.32 1s1 914. Stable 
Si 1.86 1.18 c» sii 
P 1.08 0.65 T32 Stable 
Sani 2.58 2.13 2.38 2.79 2.07 + 0.07 
Cl 3.56 3.85 3.84 3.82 + 0.06 


Lo 


jns. Up to now no other practical procedure has 
een devised for making estimates of electron 
ffinities of heavier atoms. A great deal of work 
tas been carried out on applications of different 
methods in the literature (43), and the results 
btained by several of them are in poor quanti- 
ative agreement. 


Firstly, in Table 6.1 we present results attainable 
y different methods together with available ex- 
erimental values. 

We note at once that each method gives predic- 
ons in error by a fixed constant. The constant 
> chosen so that the average error compared with 
perimental values available, is equal to zero. 
“hese constants are for the methods of quadratic 
trapolation, Bates-Moisewitch and Bacher-Goud- 
mit +0.51, — 0.42 and — 0.02 ev, respectively, 
vhereas the method of Johnson an@ Rohrlich 
orks with a remainder estimate of 10 percent. 
ortunately, advances in experimental techniques 
recent years have made reliable quantitative 
Measurements of EA(C), EA(0) and EA(S) (44). 
“he affinities EA(F) and EA(C/) have previously 
been found satisfactory by electron impacts (45). 
Some interesting points, however, are apparent 
rom data listed in the table. Thus relative affini- 
Les predicted by the method of Bates-Moisewitch 
poree within experimental errors for the 
ns F—, $— and C/—. Exceptions are the ions 
3—, C— and 0—. It is seen that B— is predicted 
instable which is contrary to the experiment. On 
ne other hand, the results obtained recently by the 
nethod of Johnson and Rohrlich (46), give 
xtremely good agreement with the experi- 
nent except for $—. Obviously, this discrepancy 
3 due to incorrect spectroscopic data which are 
wailable for the sequence iso-electronic with Cl, 
Undoubtedly, the value of EA(0) is most precisely 
jetermined by experiment, and this value also 
igrees quite closely with that obtained by the last 
nethod. [ 

The method of Bacher-Goudsmit (Hellmann- 
Mamotenko) gives results being in poor agree- 
ent with the experiments. This is also expec- 
d since the Bacher-Goudsmit theory express- 
ag an atomic state in terms of various al. 
owed energy states of the preceding ions, really 
eglects the important angular correlation effects. 


his method has been applied to many negative 
ons in the first and second period (47). The 
uadratic extrapolation procedure based on a for- 
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mula first proposed by Glockler (43), gives 
also results being quite different from the experi- 
ments. Johnson and Rohrlich (46) have 
studied the failure of the latter method and other 
extensions of it (48) to yield correct detachment 
potentials of negative ions, and found that it to a 
large extent lies in the extrapolation formulas 
themselves and a non-conscientious smoothing of 
the experimental data. The formulas have to be 
justified on physical grounds, i. e., the require- 
ment is that at least two parameters of the formulas 
must have direct physical interpretation. Evidence 
on this is given by solving the spin-less Schró- 
dinger equation for an atom or its iso-electronic 
sequence. 

By using the conventional perturbation technique 


A 
on the operator Z1 Hop — Z ! Y —, the result is 
iy Jy 


oo 
E(Z)= £ Z™ E, 
v=0 


(6.7) 


The derivation of this result is beyond the scope 
of this section but we will mention that is was 
first deduced by Hylleraas in his pioneer 
work on the ground states of the helium sequence 
(1). Replacing now E, by— a, Z? and E(Z) by 
—J(Z) as being the ionization or detachment 
potential, we have the conventional extrapolation 
formulas (49) 


JD Ea 2 (6.8) 


where a, denotes parameters to be determined. 
For the helium sequence seven of these parameters 
have been found accurate by the variational method 
(50). The labour involved in applying similar 
procedure to more complicated electron systems 
would be prohibitive. However, it is apparent from 
equation (6.8) that extrapolation formulas based 
on polynomials of higher degree than the second 
(48) or including fractional powers, have no 
physical justification at all. Obviously, semi-em- 
pirical procedures based on truncated sets of the 
infinite series (6.8) are to be preferred. In ad- 
dition, an ideal extrapolation formula must obey 
the condition of correct asymptotic behaviour when 
Z is large. This condition is 
(Z= consti E 2) 
where e is the screening constant. Imposing now 
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this condition on the series (6.8), the best choice 
of extrapolation formulas should be truncated sets 
of the series 


oo 
J(Z) =, Lif 4+-5,+ > by Lez (6.9) 
v= 


where Zepp =Z — e is the effective nuclear charge. 
It is seen that the physical constant e is uniquely 
defined by equating b,=0. The latter series (6.9) 
converges much quicker in cases where the effec- 
tive nuclear charge for the outermost electron 
Zett > 2, whereas the convergence in cases of the 
lightest ions where 1 < Zepp <2 is not so ob- 
vious. In the case of H— it is plausible that the 
series (6.9) will diverge. It is noticeable by inspec- 
tion that this will not happen for the series (6.8) 
when Z=1. Thus the method of Johnson and 
Rohrlich which is based on a truncated set of 
the seires (6.9), is not applicable to the four 
lightest negative ions. 

There remains, however, the problem of de- 
termining the parameters a, or by. This can be 
overcome numerically by using the least square 
methods, but, it is important to bear in mind that 
the first three parameters Į (45, 4, and a, in (6.8) 
and bo e and b, in (6.9)] must de determined 
from the asymptotic conditions associated with 
large Z of iso-electronic sequences. The spectro- 
scopic data of the sequences iso-electronic with the 
negative ions must be carefully checked by smooth- 


ing. 


Chapter 7. Experimental research. 


7.1 Introduction, 


It is an accomplished fact that F—, C/—, etc. 
predominate amongst the negative ions. The af- 
finities of these halogens were already experiment- 
ally determined in the early thirties by the con- 
ventional mehods based on electron impact, elec- 
tron attachment and equilibrium measurements 
(45). However, these methods failed when app- 
lied to atoms associated with smaller affinities 
and even more so for determining molecular elec- 
tron affinities. The experimental information on 
these important physical quantities is still meagre, 
and the results available are often contradictory 
for special atoms. Fortunately, the improved semi- 
empirical estimation of affinities during the last 


decade has coincided with a rapid development 
of refined techniques. It is therefore to be hoped 
that reliable experimental data on electron affini- 
ties of atoms being most important in the periodic 
table will soon be availalle. Concurrent with 
these advances there has been developed a con- 
siderable interest in the exact mode of formation 
and destruction of negative ions in different target 
gases (51). 

This chapter deals with a new experimental 
method based on the photo-detachment process 
(Section 2), and outlines principles necessary to 
understand current research on the formation and 
destruction processes of negative ions (Section 3). 


7.2 The method of photo-detachment. 


Recently, a research group at the National Bu- 
reau of Standards, Washington, has devised a new 
method of measuring atomic electron affinities by + 
observing the optical threshold for electron de- - 
tachment (52). In fact, it is safe to say that the © 
bulk of experimental works on this problem 
published during the last decade, has been done 
at this research location. The method is based on 
the fact that negative ions eject electrons by ab- 
sorption of light quanta having energy greater - 
than the electron affinity. This process may be. 
illustrated by the equation 


Z +hv+Z°+e 


(7.1) 


In principle, the problem of determining the | 
electron affinity is equivalent to finding the energy 
(ev) of the threshold of the photo-detachment - 
cross-section. As a result we have | 


EA(Z)=hc|Aq (7.2) 


where A, is the wave length of the threshold. 

The principal elements connected with a ne- 
gative ion photo-detachment apparatus are illu- 
strated in Figure 7.1. 

However, it is worth while mentioning that by 
using a specific system of optical filters Bran s- 
comb ef al. (53) obtained the cross-section to 
yield the simple expression as follows B 


o(A) =vLE(v|photon) 
XI (electron) |W (X) I(ion) (7.3) 


where E(v|photon) —hc|X is the photon energy; 
v the ion velocity, L the length of the illuminated 


“y 


y5 


2 


p=10 mm 


| Negative 
¡ion 
source 


Mass 10 
analyser 


ion 
collector 


Reaction 
chamber 


optioal 
filtors 


f 
7 


chopper 


carbon 
arc 


detached 
electron 
collector 


Figure 7.1 The principal elements of a negative ion 
photo-detachment apparatus [Branscomb et al. (52)]. 


eam, I the measured currents in the collectors 
d W(A) the radial light power. Clearly, the 
ehaviour of the cross-section in the*feighbour- 
ood of the threshold must be carefully investi- 
rated. In special cases the threshold point (À) 
an be fixed very accurately. Thus in the O— case 
he affinity is determined to an accuracy of 
1.005 ev (see Table 6.1). The assigned error 
rives confidence limits as the upper value agrees 
ith that found by Johnson and Rohrlich 
semi-empirical way. 


3 Formation and destruction of negative ions. 
Charge transfer collisions. 


It is known that the negative ion density plays 
n important róle in the inelastic collisions involved 

the ionosphere. The collisions between charged 
nd neutral particles in the ionized layers lead to 
ecombination of positive ions with electrons, to 
utual neutralization of positive and negative 
ns, to formation of negative ions by electron 
ittachment and to destruction of the same by 
‘lectron detachment. Apart from measurements 
ade by rocket-borne mass spectrometers, the ne- 
rative ions thus formed in these layers cannot be 
»bserved directly because of their complete lack 
of optical spectra. The investigations therefore 
ave to be carried out in the laboratories. 

This section is devoted to giving a short ac- 
count of principles connected with charge transfer 
ollisions. It has become customary to measure 
he charge transfer collision prohability by a cross- 
ection o(/f) where 7 and f denote the initial 
nd final charge of the beam particle, respectively. 
he cross-sections which depend on the incident 
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particle energy, can be measured separately for 
each particle in ion collectors with sensitive de- 
tectors. The fact is that each particle after charge 
changing can be immediately removed from the 
main beam by strong transverse electric or mag- 
netic fields and directed to collectors. Then the 
measured cross-section is termed as the attenuation 
cross-section by charge transfer. Clearly, the dif- 
ficulty in investigating the beam composition in- 
creases with the number of the charge components 
in question. In the current research terminology 
(54) the charge transfer collision processes are 
said to involve three charge states H^, H° and 
H+ in hydrogen and four in helium, namely 
Fez Ele Hess mand Hean 

We proceed by defining a quantity p indicating 
the number of target gas atoms per cm? which 
the beam has traversed in the pick-up tube 

p= AT K NIR T (7.4) 

where A is the Avogadro’s number, L the length 
of the pick-up tube, N the number of gas atoms 
per molecule at average pressure K dynes/cm?, R 
the gas constant and T the absolute temperature. 

However, the case of charge transfer being of 
most physical interest is that in which an equili- 
brium is established between fractions of diffe- 
rent charge components. This is the situation 
in the beam after undergoing a sufficient 
number of charge transfer collisions so that 
no further charge changing in the composi- 
tion can be detected, i. e, when p—> oo. The 
equilibrium fractions Fp(j) of ions with charge 
j cantbe computed from the observed w(i|f)s. 

In the case of Æ the fraction at equilibrium 
may be written (54) 


Fo (—1)—2ío(0|—1)o(Q-1]0) 
+ a@@(+1|—1)} 6" 
a — o (0| —1) -- o (0 |- 1) 

B= (0|+1) {w(—1]0) 
+o(-11+40)+0(+1|-1) 
+0(0|-1)(0(4+1|-1)+06(+1]0) 
+o(-1|+1)}+ 41) 0) (o (+110) 
+ (-1|+1)}+(-1|0)#(41]-1) 


For ions approximately considered as two-com- 
ponent systems, as being realized by the charge 


(7.5) 


where 


and 
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composition He and He* at kinetic energies be- 
low 80 kev [He— 0.1 percent and Het 0.2 per- 
cent] and H and H+ in energy range above 
50 kev [H— is negligible] we have the very 
simple relations as follows 


F (0) =@(+1| 0) {w(0|+ 1) 
+ w(+1]|0)* 

Po (4-1) = @(0|+1) {w (0|+1) 
+ c (+1] 0)}? 


It is seen that both attachment cross-sections 
[f «5 single and double electron capture] and 
detachment cross-sections [f 7» 7, single and 
double electron loss] enter the expressions. Simi- 
larly, the growth of different components of 
charge states can be investigated. 

The techniques usually used for transfer collision 
experiments are quite similar to those described 
in Chapter 9, Section 3. 


(7.6) 


Chapter 8. The negative hydrogen ion. 


8.1 Introduction. 


The most simple negative ion structure oc- 
curring in nature is realized through the rela- 
tively large electron affinity of the hydrogen atom. 
The negative hydrogen ion has been of physical 
interest almost since the advent of quantum me- 
chanics in conjunction with the theory of alkali 
hydrides (55). The stability of the H— ground 
state, however, was also early proved theoretically 
by using a truncated correlated trial function 
in the Rayleigh-Ritz variational scheme (56). 
Since then considerable attention has been paid 
to the H— ion as a non-deviable constituent of 
the atmosphere of the sun and similar cooler stars 
(57). The constituent accounts for the well known 
opacity. This phenomenon may be explained 
briefly as follows. In the atmosphere of the sun 
consisting mostly of neutral hydrogen atoms, there 
are conditions for creating a large number of free 
electrons by ionization of heavier atoms. The H— 
ions are then frequently formed by electron attach- 
ments. Further, the photons in the thermal radia- 
tion having energies a little greater than the 
threshold energy, are mainly absorbed by photo- 
detachment processes associated with the bound 


H-— ions. Further, negative ions are formed by 
renewed capture of the slow electrons so released. 
These two processes are the most important sour- 
ces of the opacity. 

Furthermore, the study of the H— eigenvalue 
problem is also vety interesting from a theoretical 
point of view since this most weakly bound struc- 
ture provides a severe test of the conventional 
approximation methods. 

This chapter is therefore devoted to a special 
study of this particular negative ion. 


8.2 Known facts concerning the normal ion 
structure. 


In the absence of external fields the normal 
ion is considered as a permanent structure. Even 
with an open shell approximation procedure using 
a symmetrized product of exponential one-electron 
functions, it succeeded in proving the stability of 
the ion in the ground state (58). The electron 
affinity became as good as 0.37 ev. The value of 
this quantity has been improved continually by 
applying different truncated sets of correlated 
trial functions in the Rayleigh-Ritz variational 
scheme (59) until these improvements recently 
have culminated in the excellent value of 0.75 ev 
(50). 

Since the discovery of H— as the principal 
source of the opacity (60) the H— photo-detach- 
ment cross-section is found theoretically and ex- : 
perimentally as well (61), and the theoretical | 
results are constistent with the experiments. How- - 
ever, for lack of an experimental detachment po- 
tential of H— the theoretical value, being an 
upper bound to the exact one, has been used fre- 
quently in the literature. So far the affinity has | 
not been determined from the photo-detachment | 
threshold in precisely the same way as for the 
heavier ions. The reason for this is probably the | 
difficulty encountered in fixing accurately the 
threshold point of the experimental cross-sections. 
Further, there is no urgent necessity for finding. 
the affinity experimentally in view of the ex- 
tremely good theoretical value. 

However, the last thirty years’ work on H— 
has accumulated a rather big literature, of which 
the extensive theoretical computations of the con- 
tinuous absorption coefficients by Chandra- 
sekhar (57) are undoubtedly the most impor- 
tant. It is worth while emphasizing that the radial 


^ 


d angular correlation effects are mostly equally 
istributed in the computed value (0.75 ev) 
hentioned above. The best radial limit obtained 
y superposition of purely radial configurations 
0.41 ev (12). This value leaves about 0.34 ev 
‘sociated with purely angular correlation effects. 
Finally, it is noticeable that the method of cor- 
elated eigenfunctions is superior in deriving the 
ffinity. The best value attainable by the 20 x 20 
aatrix solution when Z = 1, is by way of 
omparison 0.70 ev (see Paper I, Table 2). 


.3 Results of charge transfer collisions. 


The purpose of this section is to give experi- 
aental results of charge transfer collisions made 
hydrogen targets. We prefer to present the 
ta in graphical figures which are all reproduced 
com the excellent review by Allisorf (54). 
irstly, Figure 8.1 shows the variation of the 
-oss-sections with the kinetic energy in electron 
ilovolts [w(+]f) =0;/]. As can be seen, there are 
onounced maxima for all the curves. Cross- 
ctions associated with the negative ion, o (7| — 1) 
d œ(— 1|f), have such maxima in the energy 
gion 5—20 kev. 

Secondly, Figure 8.2 shows graphically the 
uilibrium fractions |7;4 = Fo (¿)), Fo (—1), 
7o (0) and Fo (+1) emerging from a hydrogen 
rget at various kinetic energies. It is seen that the 
tral fraction decreases nearly linearly at high 


L 


| Figure 8.1 The charge transfer cross-sections of 
components H , H^ and H* in hydrogen target 
] {Allison (54)}. 
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energies (150—1000 kev), whereas the proton 
fraction remains constant in lower energy region 
(3—10 kev) and increases successively at higher 
energies. The H-- fraction, as being the most 
attractive, shows a pronounced maximum at 10 kev. 
This implies a H— constituent of about 2 percent 
in the total beam. In other words, if neglecting 
large angle scatterings, a 1 ma 10 kev proton 
beam traversing a hydrogen target being suffi- 
ciently thick to assure charge equilibrium, will 
yield up to 2ya current of H— ions. 


Fractions Of H; H$ H* In 
A Hydrogen Beam Equilibrated 
In Hydrogen Gos 


50 100 


500 1000 


Figure 8.2. Equilibrium fractions of the hydrogen 
components in hydrogen target at various kinetic 
energies (Allison (54)). 


Finally, the Figures 8.3 and 8.4 illustrate clearly 
the growth and decay processes of the charge 
components as purely Ht and H— beams of 
15 kev traverse different target gases [ H;, N, and 
Ne] of increasing p, respectively, until charge 
equilibrium is established. It is to be noted that 
the variation of the independent variable p is 
made by varying the product in formula (7.4) 
denoted by LK as /(cm) +p (microns) at 20? C. 

Some interesting features concerning the be- 
haviour of the charge equilibrium processes can 
be drawn from the graphs. The processes caused 
by equal positive and negative current ion sources 
seem to be nearly inverse, disregarding the very 
rapid decay of the initial H— beam and the slower 
growth of final H— ions. It is also seen that the 
H— ion growth, as expected, is largest in the 
hydrogen target gas, whereas this yields the pro- 
ton growth in the neon target gas. 
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Figure 8.3. Growth of H and H9, and decay of 
H+ as a purely proton beam traverses the target gases 
Hy, Na and Ne at 20° C (Allison (54)}. 


Figure 8.4. Growth of HO and H*, and decay of 
H as a purely H beam traverses the target gases 
Ha, N3 and Ne at 20° C (Allison (54)}. 


8.4 Semi-empirical investigation of bound excited 
States, 


This section is devoted to discussion of results 
obtained by a semi-empirical investigation of 
bound excited states lying in the purely discrete 
region of the H— spectrum (62). Hereafter this 
paper will be referred to as III. 

Firstly, we consider some theoretical details. If 
we now ignore all effects due to angular corre- 
lation, the potential energy U, of the extra elec- 
tron in the field of the atomic electrons may, to 
a zero approximation, be written (40) 


mé «zo xexpi-2Z()5) (8.1) 


where Z(e) is the effective nuclear charge and re 
the nuclear distance. Obviously, this potential re- 
presents a field of force falling off exponentially 
with distance and very much more rapidly than 
a conventional Coulomb field. The angular cor- 
relation effects, of course, will in a first approxi- 
mation modify the shape of the potential curve, 
a modification being dependent individually on 
the atomic fields. 


However, from this consideration it is conven- 
ient to conclude that there may exist only a - 
finite number of bound discrete states in fields ` 
represented by negative tons. Conversely, it is well | 
known that Coulomb fields can create an infinite ` 
number of such states. Of course, natural questions 1 
as follows will arise. Does the number of discrete - 
states remain infinite when Z in the atomic Ha- i 
miltonian goes into singular points as Z — 1? 
If so, the discrete states of the negative ions will 
lie very close together and to the singly ionization 
limit. This assumption is only suggestive and needs — 
an exact mathematical investigation. So far, too 
little is known theoretically about the behaviour ' 
of the exact potentials occurring in negative ions. 
It seems plausible that there may exist negative 
ion fields creating attractive forces which can | 
make several bound discrete states possible. On 
the other hand, the current experimental technique 
is not so refined that an eventual discrete negative - 
ion spectrum can be observed, only a continuum : 
ranging from the long wave length threshold. 
Sometimes peaks occur at the threshold which 
are probably due to resonance effects arising from 
bound excited states (63). In practice we consider 
the negative ion spectra as consisting of continua. 
only and no bound energy levels in the discrete | 
region, or at most one (ground state) or two. Phy- 
sically, the metastable states are most interesting. 
Thus in the case of H— much attention has been 
paid to the (1525)'5 and (1525)3S states. Extensive 
quantum mechanical computations on the latter 
have been carried out and failed to get a bound 
state (64). 

On the other hand, existing bound states of 
H— in this discrete energy region have recently 
become interesting in conjunction with the problem: 
of elastic scattering of slow electrons by hydrogen: 
atoms (65). E 

Let us now consider the semi-empirical results 
quoted in Table 1 of Paper III. It is to be borne 
in mind that the procedure used is based on extra- 


blation formulas derived from truncated sets of 
ne series (6.8). We notice that the (1525)%S 
esult predicts an unbound state consistent with 
ae theoretical computations. Similarly, the 
325:)!$ state was found to be unbound. The 
esult concerning the latter is not listed in the 
Lble, but it implies that the state really has a fair 
aance of getting bound, at all events, much bet- 
sr than the ?$ state. 
The (1:2p)?P result, however, predicts a bound 
ate lying below the continuum limit. To be sure 
ne result has been checked by determining 
ne parameters 4, in the Johnson-Rohrlich 
theme, i.e., the coefficients a), a, and a, have 
en determined by imposing the asymptotic con- 
itions associated with large nuclear charge of 
e iso-electronic sequence (Section 6.6). . Of 
urse, we cannot place firm reliance om the re- 
It thus obtained in an empirical way; "buteit may 
iimulate attempts to investigate this particular 
antum state by conventional methods. It is 
aerefore our hope to carry out a configurational 
teraction study of the H— (152p)?P state in 
e near future. We then intend to include as 
auch as possible of configurations of the type 
ym (44-1), nq] associated with the species L=1 
d S=1, where m>g+2,n>g9+1 and 
«01:12, If now the (1s2p)*P state 
pears to be really bound in H—, it means that 
inversion of the (152s)*5 and (1s2p)*P levels 
s taken place. Assuming that the (1525)!S state 
also really bound in H—, an investigation of the 
32p)!P state might be of interest to see if there 
an inversion of these levels in violation of 
nventional rules in the LS coupling scheme. It 
ould be noticed that the result obtained for the 
s2p)*P state when using similar semi-empirical 
nethods implies such a possibility. On the whole 
>mi-empirical estimations of detachment potentials 
the H— %*%P terms provide promising aspects 
getting bound states although the results may 
considered very uncertain in view of probable 
Irge errors in the observed data of the higher 
citation potentials. 


5 Quantum mechanical investigation of bound 
= quasi discrete states. 

| This section is concerned with results obtained 
y quantum mechanical computations of some 
1asi discrete energy states of H—. Most attention 
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has been focused on the lowest exactly quantized 
(2p?)*P state (66). The computation is quite 
similar to the analogous helium state (Part 2, 
Chapter 5). It is to be noticed that the same 
truncated basis is used, and here we confine our- 
selves to discussing the results only. As can be 
seen from (66, Table 1) the best values attainable 
from the optimum 10x10 matrix when Z=1, 
are Nopp = 0.3180 and Emin= — 0.248543822 
Rydbergs. It is also noticeable by inspection that 
the optimum point on the E(7) curve is determined 
very punctiliously. Thus the check by equations (6) 
of Paper I gives errors smaller than 4 parts 
in 10° for y and 7 parts in 10° for E. If we now 
compare this computed H— energy with the exact 
value of the H (2p)?P state (the ionization limit of 
—0.25 Rydbergs), we get a negative electron affi- 
nity — 0.020 ev. This value represents absolutely a 
lower limit of the detachment potential. However, 
our H— computation should have a considerable 
remainder error, probably of the order 0.030 ev. 
To make an estimation of a reasonable upper limit 
of the potential, we point out that the remainder 
error in H— is certainly not greater than that of 
the analogous He case. Let us now suppose this 
error to be of the same amount. As previously 
mentioned in Part 2 (Section 5.8), we may 
in fact consider two experimental values of the 
He (2p?) °P state. In this way we deduce total 


H— energies of —0.2503 Rydbergs (£46. = 


—1.4220 Rydbergs) and —0.2573 Rydbergs 
(EDS, =71.4290 Rydbergs). These values 


correspond to electron affinities of 0.0041 and 
0.012 ev, respectively. Obviously, we may consider 
the latter as representing an absolute upper limit 
of the detachment potential. From these considera- 
tions it seems apparent that the (2p?)*P state is 
bound in H—. The state is not metastable as pre- 
viously assumed in the communication (66) since 
it will decay by allowed dipole transitions to 
probably lower bound H—  (1szp)?P states 
(Paper III). 


Until now we have neglected completely spin- 
orbit and spin-spin interaction effects in our 
treatment of the exactly quantized quasi discrete 
states. As emphasized previously, corrections due 
to such effects must be taken into account for 
treating heavier atomic systems. Let us see how 
this works in the case of the two-electron atomic 
(2p?) *Po, 1, 2 state. It is clear that through partial 
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breakdown of the LS coupling in the (2p)* con- 
figuration, the two components 3P, and *P, become 
mixed with the imperfectly quantized states *So 
and 1D,, respectively. The third component *P, 
is not subject to such breakdown and remains 
therefore exactly quantized. The transition prob- 
ability of the *P, component should be of the 
order 10% sec— (only dipole transitions to lower 
states). On the other hand the other two compo- 
nents thus being subject to autoionization, have 
transition probabilities of the order 1013—1014 
sect. At all events it is seen that a bound 
H— (2p?)*P state should have a lifetime long 
enough to account for observation, but, on account 
of the very small electron affinity, this seems to 
be unfeasible. 

Similar theoretical computations have im 
three imperfectly quantized states being bound in 
H—, namely the states (25?)'S (see Paper I, 
Tables 1 and 3), (2:2p)?P (66) and (3p?)*P 
(66). The most striking feature about the results 
obtained is that the detachment potential in the 
(252) 18 state (0.82 ev) is larger than the ground 
state electron affinity (0.75 ev). 

A common feature of the imperfectly quantized 
states is that they, roughly speaking, have a larger 
tendency to get bound. Apparently this effect is 
due to the strong interaction with the adjacent 
continuum. Physically such bound states in negative 
ions are less important because of their short life- 
times (10—%—10—" sec). It is worth while men- 
tioning that they are interesting in conjunction 
with resonance states of different compound 
ions (67). 


8.6 Conclusion. 


From the preceding discussions it is obvious 
that the problem of the most simple negative ion 
structure is still far from being solved completely. 
It is true, of course, that its existence in the ground 
state has been proved, but the theoretical search 
on bound excited states having discrete energies 
in the H— spectrum is still in its infancy. 

We have seen that the (1:2p)?P state seems 
highly probable to be bound in H—. Work on this 
problem is being planned, and probably a positive 
result may one day become significant for develop- 
ing a more reliable theory of the H— ion. Further, 
our computation indicates the possibility of a 


bound (2p?)*P state which thus becomes ex- 
tremely good as initial upper level for dipole 
transitions. A weakness of the theoretical search 
up to now has been that too much effort has been 
concentrated on purely atomic and molecular pro- 
blems, whereas negative atomic ion systems have 
been scarcely treated. It is to be hoped that a com- 
prehensive extension of our (2p?)*P computation 
will be carried out in the near future. It is likely 
that improved approaches may one day change a 
little the framework of the present view on the 
H— structure. 


Chapter 9. The negative helium ion. 


9.1 Introduction. 


In this chapter we deal solely with the negative 
helium ion which recently has been of renewed 
interest. 

The presence of a minute constituent of He— 
in canal rays from a purely helium discharge was 
observed by H i by in 1939, as a particular peak 
in the mass analyzed spectrum (68). In fact, the - 
earliest evidence of a He— constituent dates back 
to a published paper by Dópler in 1925 in- | 
volving a similar analysis of canal rays ejected 
from a mixed discharge of hydrogen, helium and 
neon (69). Unfortunately, no reliance was placed 
on this result since the He— track on the screen 
was assumed to originate from impurities or oc- 
currence of some recombination processes in the: 
target gas mixture. 

From a theoretical point of view, however, even 
many years after Hiby's detection, it was believed 
that He—, one of the least stable electron-negative 
ions, could not exist as a permanent structure. 
Nor wete further experiments in producing the 
ion made. Apparently, the detection of a probable 
He— ion stimulated several theoretical attempts 
to find any singly excited state being bound in 
the ion (70), but without result. It was easily 
proved that the ground state (1:22:)?$ was un- 
bound. Fortunately, this futile search on bound 
singly excited states resulted in focusing con- 
siderable attention on the quasi discrete states of 
the ion (71). To explain the existence of He— 
in high vacuum under special experimental con: 


| ditions, it was early suggested that a possible 
bound quasi discrete state must be metastable in the 

¡sense of the He (1525)3S state. Therefore atten- 
tion was concentrated on the (15252p)4P state 
being the lowest quartet. Extensive quantum 
mechanical computation of the total He— energy 
in this state indicated a positive electron affinity 
of 0.075 ev (72), i. e, the state is obviously 
bound. This paper published in 1955 will here- 

| after be referred to as IV. 


Since then the question of the He— existence 
has been approached by a number of experiments 
which, in fact, have added greatly to our meagre 
knowledge of properties of this peculiar negative 
ion. 


9.2 Experimental evidence for the extiience 
of tbe ion. 

Compared to H—, so far little is known about 
properties of He— obtained by experiments. In 
the period 1940—1955 no attempts were made 
to be quite sure of its existence or not. It is 
highly probable that the reason for this is to be 
found in the prevailing view on the He~ non- 
‘existence which still was in evidence among the 
majority of the experimental physicists. Ap- 
parently, the theoretical proof of a bound meta- 
stable state in He— (Paper IV) and the develop- 
iment of vacuum techniques in the intervening 
¡years have stimulated research groups to verify 
. Hiby's detection (68). 

In fact, the investigations start with a short 
report of Dukelskii eż al. (74) which quotes 
observed cross-sections of He— in different target 
¿gases. Immediately on this publication Win d- 
ham ef al. reported quite briefly a pronounced 
e— peak in the observed momentum spectrum 
of different negative ions (73), apparently quite 
ignorant of the Russian report. In a subsequent 
aper Windham ef al. (75) have given a 
etailed description of the He— investigation. 
The apparatus used is schematically shown in 
igure 9.1. 

A detailed description of the technique used 
ere is beyond the scope of this section. We only 
motice that in the analyzer the various components 
are spread out by magnetic fields into a momentum 
spectrum showing significant peaks of present 
egative ions. 


ps 
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We reproduce two He— graphs from the work 
reported by the American group (75). Figure 9.2 
shows the momentum spectrum of negative ions 
when the incident particles are 17.5 kev Het ions, 
and the target gas is hydrogen. As can be seen, 
the spectrum provides conclusive evidence of the 
presence of He— ions. Figure 9.3 shows the He— 
current as a function of the extractor potential. 
It is seen that the current increases with the 
voltage. 

The group succeeded under these experimen- 
tal conditions in getting the He— ion current in 
amounts as large as 3.5 x 10—? amp. This value 
has been improved by another American research 
group to 5 x 10^ amp when using a more refined 
He* source (76). 
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Figure 9.1. Schematic diagram of the apparatus. 
{Windham et al. (75)}. 
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Figure 9.2. The negative ion momentum spectrum 
with He* ion source and hydrogen target 
{Windham et al. (75)}. 
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Figure 9.3. He current as function of the extractor 
voltage {Windham et al. (75)}. 


9.3 Charge transfer collision experiments. 


Recent years have also been characterized by 
experimental activity in helium charge transfer 
collisions (54). In fact, we may consider helium 
as a four component system consisting of the charge 
states He—, He", Het and He?*. As usual, the 
ion currents are measured by special detectors 
placed in separated ion collectors. The cross-sec- 
tions and the equilibrium fractions can then be 
computed as previously explained in Section 8.3. 

The Russian group (74) has determined at- 
tachment cross-section o(--1|—1) arising from 
two-electron capture in target gases as helium, 
neon, argon and krypton. This quantity is plotted 
graphically in Figure 9.4 in the kinetic energy 
region 15—175 kev. It is to be noted that the 
units used are 10—!? cm? per target atom. Ob- 
viously, the cross-section depends appreciably on 
the target gas being used. It is seen that for argon 
and krypton targets the picture seems to be quite 
similar as to the shape of the curves and with 
maxima at nearly the same place (65 kev). For 
the neon target the maximum occurs at 70 kev, 
whereas no maximum is observed in the helium 
target because o(4-1|—1) seems to keep constant 
independent of the incident Het beam energy. 
This result seemed at first to be a little surpris- 
ing, but the group concluded that a determination 
of the exact shape of the He curve was impossible 
with the apparatus used, since the cross-section 
in that case came out with a. maximal order of 
10—21 cm?. The shape, however, is accurately de- 
termined by an American research group (77). 
At all events, the results obtained by the Russian 
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group are important as they indicate clearly that 
the He— ions are frequently formed by double 
electron attachment in single He* collisions with 
atoms of the noble gases neon, argon and krypten. 
This assumption has been supported by later in- 
vestigations (77). The growth of He— in aa oti- 
ginally pure Het beam as it traverses the hydro- 
gen target, is graphically illustrated in Figure 9.5. 
The solid curve is drawn to give the best fit with 
the observed data. It is seen that the growth up 
to the maximum is nearly linear indicating really 
that o(+1|—1)>>o(0|—1)>0. Comparing the 
result with that of the analogous case of H— (see 
Fig. 8.3) we note a striking discrepancy. We see 
that the growth of the He— ions first reaches a 
broad maximum and then decreases with increas- 
ing gas flow in the pick-up tube. Conversely, the 
H— growth reaches saturation with increasing 
gas pressure. Clearly, this discrepancy may be ex- 
plained by the different behaviour of the consti- 
tuents of the negative ion beams. In the H— 
beam the ions occur frequently in the stable 
ground state, whereas in the He— beam the ions 
are formed in any metastable state limited by life 
time of finite order. This conclusion is further 
supported by extended investigation of the He 
growth in other target gases (77). The equilibrium 
data of He—, however, are accurately determined 
in a few cases of target gases (77). The determi- 
nation of the helium equilibrium fractions 
F (+1) and Fe (—1) at various He* energies 
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Figure 9.4. The cross-section wH | —1) for Bea k 
moving in target gases He, Ne, Ar, and Kr. i 
{Dukelskii et al. (74)}. 
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has been made recently by Jorgensen et al. 
(78). A He* beam is sent through a hydrogen 
target at constant 30 » pressure and of variable 


| length to ensure charge equilibrium in each case. 


The equilibriated beam emerging from the pick-up 
tube is analyzed in a vacuum magnetic field. Thus 
the Fo (—1) fraction is determined in the kinetic 
energy region 40— 200 kev. Corrected data ob- 
tained by the same group (77) when using im- 
proved experiments are listed in Table 9.1. 


Table 9.1. Equilibrium data concerning the 
He ion. Maximal values. 
, Kinetic 10) € Kinetic 
Target gas BO energy — Fs (—1)/ energy 
(kev) FR (+1) (kev) 
H, 2.3 130 ME 05 
He 0.95 200 25 165 
N, 7, 80 DED 60 
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Figure 9.5. The growth of He in an originally pure 


He* beam as it traverses the target hydrogen. 
{Windham et al. (75)}. 


9.4 Lower limit of the lifetime of the ion. 
We explained the peculiar behaviour of the 


| He— growth in charge transfer collisions as being 
| due to ions excited in any bound metastable state. 
| As is well known, if the investigations are not car- 
| fried out in absolute vacuum, the lifetime of a 
| metastable state is by no means infinitely long 


but limited by collisions of the second kind with 
the residual gas atoms. In general, usual dipole 
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radiation from such a state is quite impossible. 
A metastable He— ion will thus decay solely 
through electron collisions or atomic collisions of 
the above mentioned type. Its high excitation 
energy will then be used in exciting residual gas 
atoms. Theoretically an isolated He— ion should 
have an infinitely long lifetime, but even if we 
could produce a pure He— beam in isolation tech- 
nologically by using a very high vacuum drift 
tube, a natural He— decay would occur because 
of the existing small detachment potentials. It is 
plausible to think that by continuous inter-He— 
collisions quite a number of slow electrons will 
be released from the ions. A dissociation process 
as follows may then be possible 


He estou n He” + 2e (9.1) 


From this consideration it is obvious that even a 
He— beam in isolation must have a lifetime of 
finite order. However, it has been observed that 
metastable constituents have a lifetime up to the 
order of 1 sec. On the other hand, it is clear that 
the lifetime of He— in gas discharges is large 
compared to that of normal excited states of about 
the order of 10—* sec. Of special interest is, of 
course, the lifetime arising from natural decay. 
To determine this exactly from experiments using 
the conventional vacuum technique, is not possible, 
but a lower limit can be found approximately. 
This has been done quite recently by Sw eet- 
man (79). The apparatus used consists mainly- 
of a 60 cm long drift tube, schematically shown 
in Figure 9.6, into which a pure He— beam is 
injected. At the other end of the tube the result- 
ing components He? and He- are analyzed by 
magnetic separation. The lifetime associated with 
natural decay may then be determined from the 
simple relation 


t=L]v In{1 — (1--a[B)—] (9.2) 


where L is the length of the drift tube and v the 
particle velocity. It is to be noted that œ is the 
number of He? particles produced by natural decay 
of the He— beam along the 60 cm long flight 
path in the tube, and 8 is the number of particles 
in the analyzed He— beam. The pressure main- 
tained is less than 10—* mm Hg. A lower limit 
of r is found as follows. The He” component 
measured at the end of the tube arises essentially 
from natural decay of He— together with dis- 
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sociation by collisions with the residual gas. Dis- 
sociation effects as being due to diffusion to dif- 
ferent parts of the wall and beam stop are estimated 
to be negligible. If now equating « in (9.2) 
with the total observed He? component, we get 
an upper limit of the probability of natural He— 
decay and thus a lower limit of 7. The minimum 
lifetime of He~ found in this way is 1.0 x 107 
sec. In the light of this result it is easy to under- 
stand why recent attempts to obtain high power 
He— ion sources were doomed to failure when 
accelerating pure He~ beams in Van de Graff 
generators with conventional vacuums. The loss 
of He— ions along the flight path will be too 
large. One has therefore to operate with much 
higher vacuums. It is thus found that the residual 
gas pressure must be reduced to below 5 x 107? 
mm Hg for keeping the He~ loss less than 20 
percent (79). 

Of special interest is a comparison with the 
lower limit obtained experimentally for the me- 
tastable He ?$ state (80). This limit is of the 
order 103 sec, and is thus much larger than that 
of He—. Since the detachment potential of the 
helium state is about 45 times higher than that 
of He— (Paper IV), this fact seems to be reason- 
able. 
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Figure 9.6. The He drift tube {Sweetman (79)}. 


9.5 The existence of the ion in the metastable 
(15252p)*P state. 


Extensive quantum mechanical computation of 
the. total He— energy in the metastable 
(1:2:2p)*P state has been carried out. Compu- 
tational details are given in Paper IV. In this 


section we shall discuss the final result which 
predicts a bound He— state, namely 
Esop (He 35) — EHe *P) —0.075 ev (9.3) 

Clearly this value represents an absolute mini- 
mum of the detachment potential since the com- 
putation is based on the Rayligh-Ritz variational 
method. It is highly probable that the computed 
He— energy is still in considerable error since 
only angular correlation effects as being due to 

rj —terms in the eigenfunction have been taken 
into account. This error is certainly not greater 
than E,(He?$)— Eegpti (He *S)=0.075 ev (see 
Paper IV) where E, is obtained by a two-electron 
eigenfunction of the same form and flexibility. 
This gives a reasonable upper limit of 0.150 ev 
for the detachment potential. In fact, we feel that 
the average value of 0.11 ev is quite close to the 
exact one. It is thus seen that the normal H— ion 
should have about a seven times higher detach- 
ment potential. This fact certainly accounts for 
the different behaviour of the H— and He— ion 
in charge transfer collisions. 


The *P state lies in the quasi stationary energy 
interval of the He— spectrum. It is easy to see 
that the state is exactly quantized when using 
the conventional LS coupling scheme. If, making 
allowance for spin-orbit interaction, we see that 
the component ‘P,), is still exactly quantized, 
whereas the components *P,, and *P,, being 
capable of a breakdown of the LS coupling, be- 
come mixed with the imperfectly quantized states 
?p,j, and ?Py, associated with the same con- 
figuration, respectively. 

Furthermore, the metastability of the (15252p)*P 
state is obvious since it is the lowest quartet of 
the ion. No dipole radiation can take place from 
this state. The decay probability of the 4P,), 
component will thus depend on collisions of the 
second kind with neighbouring particles. In iso- 
lation, however, such effects should be negligible. 
In fact, if making allowance for magnetic spin- 
spin interaction between various configurations 
in a second-order perturbation, the *P,,, level is 


also subject to autoionization. The radiation-less — 


transition (15252p) *Py, —(1:2&f) 27, can take 
place consistent with total angular momentum 


and parity conservation. So far it has been thought — 
that the spin-spin operator possesses negligible - 
elements for this transition in He—, but quite 


v 
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recently (96) the autoionization probability as- 
sociated with this, has been calculated to be of 


| the order 10? sec—!. As can be seen, the result 


agrees very well with the experimental value (79) 
as indicating a lower bound of 10—* sec for the 
He— lifetime. In conclusion, it is worth while 
mentioning that the imperfectly quantized state 
(2:2p?)*P is found by similar quantum mechani- 
cal methods to be bound in He—. However, the 
result is somewhat uncertain as the difference 
E.(He— *P)—E,(He—*P) is found to be po- 
sitive (71). 


9.6 Results of semi-empirical investigation of 
bound singly excited states. Discussion. 


Semi-empirical search on boufid states in the 
purely discrete region of the He— spectfum has 
also been carried out (Paper III). The extra- 
polations are based on the same formulas as used 
in the H— case (Section 8.4). As can be seen 
from condensed He— data i Table 1 of Paper III, 
the results predict an unbound ground state con- 
sistent with performed quantum mechanical com- 
putations (81). Furthermore, they predict bound 
He states, such as (15 25)'P and (15° 3p)?P. 
The result associated with the latter, however, is 
very uncertain in view of probably much larger 
errors connected with observed data of the higher 
excitation potentials. It is very likely that the 
former is actually bound in He—. If so, the com- 
puted total He— energy of this state must be 
E,(He— 2p °P) <— 5.8074 Rydbergs. At first 
sight such a low expectation value of the energy 
seems to be unattainable by using any correlated 
eigenfunction, but it may be possible from the 


following consideration of the problem. 


Let us as a first approximation neglet the 
mutual interaction of the K-electrons. Contribu- 
tions to the detachment potential then originate 
from probable binding energies in the purely H— 


states (152p)*P and (1:2p)?P. It is to be borne 


in mind that both these states are probably bound 
in the H— structure. Whether these component 
states can result in binding the He— structure, will 
to a high degree depend on the mutual inter- 
action of the K-electrons. If now the 2p?P state 


| is bound, it is at all events clear that this state 


cannot account for the observation of He~. We 
have seen that the observation requires a minimum 
lifetime of the negative constituent to be of the 
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order 10—5 sec. Obviously, dipole radiation is 
possible from the 2p?P state, and allowed transi- 
tion as 2p?P — 2;?$ can occur. Thus. in that case 
the He— ion will decay after about 10—* sec (the 
ground state 25?S is unbound). 


9.7 Conclusion. 


While an attempt has been made to give a 
coherent picture of the present state of the He— 
research, the main emphasis has been placed on 
problems arising in conjunction with the existence 
of the ion by experiments. We have seen that the 
ion from a quantum mechanical point of view 
can exist in a mestastable state with a lifetime 
long enough to be detected experimentally. This 
is in contrast to the previous conventional belief 
in its theoretical non-existence. The detachment 
potential of the ion being so far estimated to 
0.11 ev, is about one-seventh of that of the nor- 
mal H— ion. It is clear that this relatively small 
potential effects in shortening the He— lifetime. 
Apparently, the He— ion will represent a minute 
but observable constituent in pure helium beams 
as long as the vacuum technique is not completely 
developed. If such technological capabilities could 
be possible in the future, high power He— ion 
sources will certainly be very useful. The problem 
of existing He— ions in the atmosphere of dif- 
ferent stars and in planetary nebulae of extremely 
small density, might be very interesting. At all 
events, the He— reseatch is still in its infancy. 


Chapter 10. The negative lithium ion. 


10.1 Introduction. 

Negative ions formed by the most electro- 
positive elements such as hydrogen and the alkali 
metals L;, Na and K ate of exceptional interest 
and have been studied by a number of investi- 
gators. In Chapter 8 we have given a brief survey 
of the current situation of the H— research. We 
now proceed to outline quite briefly the present 
state of the L;— research. 

Firstly, it is worth while mentioning that 
pioneering works on the. existence of the Na— 
and K— ion carried out in the thirties (82), used 
numerical SCF functions including exchange, and 
the results were negative. This failure may be 
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ascribed to the fact that angular correlation effects 
are completely neglected in the conventional SCF 
scheme. In the meantime the existence of these 
ions has been confirmed by discharge experiments 
(83). 

The negative lithium ion was first observed by 
Sloane and Love (84). It was apparent from 
the very beginning that the observation of L7— 
was prevented by incapabilities of the vacuum 
technique used. On the other hand, the theoretical 
computations were connected with special mathe- 
matical difficulties when treating this four-elec- 
tron atomic system. With the development of the 
modern electronic computers in the future it is 
to be hoped that a thorough quantitative determi- 
nation of the electron affinity is possible. So far, 
this has not been accurately found. Different 
values ranging from 0.38 to 0.82 ev are suggested 
(85). Curiously, the average value of the above 
mentioned limits (0.60 ev) seems to be close to 
the exact one (86). 


10.2 Experimental research. 


Comparatively little is known about the L7— 
ion by experiments. Thus very little experimental 
data have been obtained on charge transfer col- 
lisions involving Li— ions (94). The maximum 
fraction /.(—1) in the target gas N, is found 
to be 2.3 X 10^* at about 150 kev. 

As shown in Section 7.2, it has become possible 
to measure the photo-detachment cross-section in 
the laboratory and thus make quantitative deter- 
mination of the electron affinity by fixing the 
threshold point. As previously mentioned, such 
experiments have been made for H—, and the 
results obtained were in fair agreement with those 
found theoretically (61). On the other hand, as 
to our knowledge so far the L7— photo-detachment 
cross-section has not yet been found experiment- 
ally, only theoretically when using the empirical 
value of 0.38 ev for the L7— detachment potential 
(87). Apparently the reason why the electron 
affinity of neither H nor Li has not been deter- 
mined with equal accuracy as that of C, O and S 
(44), is that special difficulties were encountered 
in observing exactly the optical threshold for the 
electron detachment. In the absence of such an 
experimental L;— detachment potential, we are 
in fact forced to compute it theoretically either 
by pure quantum mechanical or semi-empirical 
methods, or by any combination of both. 
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10.3 Semi-empirical research. 


Empirical estimated electron affinity of 0.38 ev 
for atomic Li is found by direct quadratic extra- 
polation of spectroscopic data of the Be sequence 
isoelectronic with L;— (45). This value has been 
improved to 0.39 ev when using a higher order 
of the extrapolation (48). The latter extrapola- 
tion, however, has been criticized by Johnson 
and Rohrlich (46). As previously emphasized, 
the extrapolation formulas must be justified on 
physical grounds, and errors occurring in observed 
data must be corrected by smoothing before the 
extrapolation. We have done this in the Lz— case, 
and the empirical value is improved to 0.476 ev 
(88). The result is consistent with a combined 
quantum mechanical and semi-empirical computa- 
tion (89). As we shall see later, this value also 
turns out to be a lower limit of the probably 
exact one. 

Of special interest is also any singly excited 
state which may be bound in L;—. A similar 
estimation has therefore been made empirically 
on the states (1522p2)1D,  (1522p2)S and 
(152252p)**P (Paper III), and the result predicts 
a bound 2D state. It is of interest that the maximal 
value of +0.367 ev for the detachment potential 
is deduced by using the extrapolation formula 


Iy(Q=a277+a,Z2+4, 


+a, Z +a, Z (10.1) 


10.4 Quantum mechanical research. 


Several quantum mechanical attempts partly in 
conjunction with empirical estimations have been 
made to compute the detachment potential of Lz— 
or the electron affinity of L; (90), but with re- 
sults not being very satisfactory. The shortcomings 
of these attempts are due to the fact that they 
essentially neglect the very important angular cor- 
relation effects which arise from the 2s-electrons 
of different spin in the L-shell. If we in a first 
approximation ignore the  inter-shell inter- 
action in Li~ and consider the K-shell as being 
mostly at the nucleus, the binding energy arising 
practically from the polarization effects of the 
L-shell, amounts to 0.81 ev (91). 

The attempt which first really proved a positive 
electron affinity of Li without using any empirical 
estimation, was made by the author in 1954 (92). - 


~~ 


e introduced a suitable chosen correlation factor 
14-«r2,]into the eigenfunction and obtained 
e difference E,(Lz *$) —E,(Li— 18) = +0.18 ev. 
.n inherent difficulty in treating the L;— ground 
rate by using correlated trial functions in the 
-ayleigh-Ritz varational scheme is that the in- 
roduction of explicitly inter-electron distances 
omplicates the basis considerably. A more prac- 
cal approach is the method of configurational 
fateraction using an entirely discrete basis, such 
8 the orthonormal standard orbital set 


WX exp (— Er) L2, 41 (251) Yon (9, p). 

The flexibility of the basis thus being reduced, 
- may be offset by the very great simplification 
evaluating the energy integrals. Up to now we 
ave only planned such a computatiópn-on the L;— 
ergy in the ground state. However, interest has 
iminished considerably since we quité” recently 
earned that similar extensive computations are in 
jrogress (86). Nevertheless, we will here provide 
Dome device whereby to attack the problem. In 
eating a four-electron system the danger lies 
ssentially in the accumulation of errors due to 
me approximation. To keep this to a reasonable 
ninimum we intend to use only one inter-shell 
Dlitting parameter (non-linear) in addition to an 
rver-all scale parameter. The energy integrals may 
expanded as linear combinations of previously 
aluated two-electron integrals (Paper I) and 
Dmputed numerically on the electronic computer 
ing similar programmes. It should be mentioned 
bere that a research group in Japan has made 
ecent configurational interaction studies on the 
;— ground state using radial open shell con- 
,gurations solely (90). As expected, the results 
came negative. On the other hand, a Russian 
xoup (89) using an extended many-configura- 
ional approximation in the conventional HF 
heme, succeeded in getting an affinity of 
.462 ev when taking into account amending 
pnfigurations such as {(15)”(2p)?} and {(2s)? (2p)*}. 
| thorough configurational study on L7 in its 
mound state using exponential orbitals of the 
mnventional Slater type NX exp [— Sr) pr 
Co KO, p) in the basis, is in progress at the 
Tniversity of Chicago (86). The preliminary result 
tained by a basis including 55 different con- 
igurations, predicts exp 4 (Li?S) — Z. (Li *S)= 
4773 ev. 2 ! 

This value is close to our empirical estimated 
alue of 0.476 ev. Since the variational method 
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is used, the value 0.4773 ev represents a lower 
limit of the exact one. By equating the remainder 
error in the L;— energy thus computed with that 
corresponding to the Be computations, Weiss 
et al. (86) deduce an upper limit of 0.6547 ev. 
If extrapolating the errors along the Be sequence 
to yield L;-, the group obtains a detachment 
potential of 0.6158 ev. The results are published 
elsewhere. 

In conclusion, from these considerations it 
appears highly probable that the value of 0.62 ev 
gives a very good approximation to the exact one. 


10.5 Conclusion. 


In view of the meagre knowledge of exact 
properties of the L;— ion, it appears as if an 
experimental determination of the detachment po- 
tential would be desirable. If we accept the theore- 
tical value of 0.62 ev, it suggests a structure being 
looser bound than that of the H— ion. In the 
first place, a careful measurement of the L;— 
photo-detachment cross-section would be of extre- 
mely good use in order to fix the threshold point 
of electron detachment as accurately as possible. It 
is to be hoped that future experimental work on 
Lr- will be centred around this problem. In the 
meantime, from a theoretical point of view, we are 
looking forward to seeing the final results from the 
American research group (86). Beyond this pro- 
blem there still remain interesting questions such as 
whether bound states exist, if at all, either in the 
purely discrete energy region or in the quasi sta- 
tionary energy interval of the L7— spectrum. These 
problems will be left to be decided by future quan- 
titative quantum mechanical computations. 


Appendix A 


The evaluation of the integrals an(a|q). 


Evaluating the overlapping integrals in (4.5), 
we restrict this evaluation to involve discrete hy- 
drogen-like orbitals only, and the case g=0, i. e., 
we consider merely the as[»s] overlappings. 

Fortunately, this may be easily done by making 
use of the contour integral representations (28) 
of the normalized orbitals 
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p*(n0| Zr/n)= — (2104)? 


n? (2n)? 
Zü d» 

X J dss” (1—35) ^ X exp s a = (1) 

DN (ao | nr) = (2m) è? {a (a + 1) (2700) 

x [drr*(1—27* X exp E sd (2) 


1» 


The contour C may here be reduced to circles 
around the poles s=0 and /—0, and we ignore 
the minus sign in (1) which is due to the con- 
ventional definition of the Laguerre function. 
Applying now the well known Cauchy's theorem 
of residues, we finally deduce 


.. 8(20 (28 R 
a, (a | 0) — ia(a-F Dj a> ES 
Xx {1 —87}” Fa (v) Gr (v) 
where Y 


Fa(v) = es (—1)e*»*1 §a-(@+1) 


G, (vy) = 1621) + (22)] §7z-0+1) fete (142) 8 n—(v+2) 
Jey 
wine 


two integers @—1 and 7z— 1. 


and o is the smaller of the 


E= Zin, ð = 


Appendix B 


1. The evaluation of M(1|j). 


When evaluating the matrix elements of the 
kinetic energy operator (5.6), we have used suc- 
cessfully the following recurrence relations 


0 1 1 2 
ap Pe la P, — Pj] 


1 
and cog Pi=> {gg +1) P +P) — (1) 

It should be mentioned that these relations are 
derived from the general and more known recur- 


rence formulas connecting the associated Legendre 
polynomials 
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2 y = Y cotg Rol and cotg DP Ld 
v, pv+2 
(2 DI! or) -0»8ETPE) 0 


Analogous with the *$ computations in Paper I 
the matrix elements are obtained with a non- 
normalized total eigenfunction of the form (1.13) 
rather than of the form (5.3). We have only to 
replace the angular part by F(9, 9, q|1, 7) given 
in (5.2). The coefficients c(m, n, q|n) in (5.3) 
are then derived from the results obtained by the 
electronic computer by using the equations (I, 14). 

Substituting now the eigenfunction whose an- 
gular part is subjected to the normalization con- 
dition f 7(9,9,,p|1,p) 7(8,0,,0|1,7) dU — 8p, into 
the expression (5.6), and making use of the re- 
currence formulas (1) in reducing the angular 
combinations and of the contour integral repre- 
sentations of the radial orbitals (I, 4), we finally 


deduce M (ab | pg | ed) — 1 (A (ab | g) 8 ca | 


+ A (cd | g) da + A (ad | q) 85, + À (bc | q) 9,4) ô pg 
where the integrals A(a@B|q) are equivalent to 
those found for the 15 configurations in Appen- 
dix (equ. 16) of Paper I. 


nmm tt 


2. The derivation of C4,Cp|&|q). 
The problem is now to evaluate the integral 


W (p|&|4) — J P} (cos 9) P; (cos 9) P; (cos 0) 
X sin 0 d9 = Cas (p | k | g) 95.552, 2s (3) 


where Ps (cos 9) comes from the expansion over 
VU 
Fortunately, Gaunt has deduced previously. 
in his work on the helium triplets (95) an explicit 
expression of the integral over the product of 
three different associated Legendre polynomials. 
Thus the expression of C,,(p|£|g) in (5.9) is 
derived from his formula (95, equ. 9) when 
necessary angular normalization conditions ha 
been taken into account. The range of y in (5.9) 
is sufficiently indicated by the factorials in the 
denominator. Furthermore, the integral W (p|k|q) 
is limited by the triangular conditions. In other 


~~» 


vords, the integral will vanish unless the integers 
|, q and & are equal to the sides of a triangle of 
[ven perimeter, i. e. 


8 Bel Peg 2s 
btgt£,2s — 10, bt gt £ZE2s 


5=0,1,2... (4) 
It is thus seen that the number & is limited by 
atisfying the condition 


pte&»|p—4| (5) 


- The matrix elements L'(1|j). 
'The sequence of interacting configurations are: 


PY”, 2p5p, (3p)”, (34), 2p4p, Sp4p, (4p), 
4d, (44)? and (4f)?. Similar to Paper T, we 
ly give numerically the elements in thé upper 
riangle of the symmetrical 10 x 10 matrix as fol- 
»ws. Each row starts at the diagonal with symbol 
noted by |mg z4| and continues to the right 
ntil all the configurations according to the above 
rder of interaction are included. 


|(2p)?| 21/16, 21/42V5, 27/160, 33/64V 5, 
5/64V 15, 81/640V3, 17/320, 99/128V 35, 
5/224V 5, 143/640V14,  |2p3p| 21/32, 
3/320V 5, — 33/640, 237/640V 3, 
59/320V 15, 111/640V5, 33/160V7, 
3/4480, — 429/1280V 70, |(35)?| 207/200, 
1/80V 5, 111/160V 15, 1743/3200V 5, 
9/3200, — 121/640V 35, 69/280V 5, 
131/6400 14, |(34)?| 9139/8640, — — 
-11/160V 3, 11/320V 15, 167/3840V 5, 
139/17280V 7, 16067/120960, *: 
63/1536V 70, |2p4p| 173/320, 11/32V 5, 
141280V15, —121/640V21, 9/4480V3, 
- 273/2560V 210, |3p 4p| 1841/3200, 
313/12800V 3, 77/320V 105, 

- 1399/8960V 15, — 1131/25600V 42, 

4p)? | 11193/12800, 361/7680V 35, 
659/26880V 5, 1413/25600V 14, |34 44 | 
1139/17280, 12431/26880V 7, 

- 3263/21504 10, |(42)*| 80543/94080, 
699/1075270, | (4f)?| 24769128672. 
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